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Abstract 

We study the motion of an inextensible string (a whip) fixed at one point in the 
absence of gravity, satisfying the equations 

Vtt = Ss{crVs), CTss - IVssl'^cr = -IVstl"^, hsP = 1 

with boundary conditions r]{t, 1) = and a{t, 0) = 0. We prove local existence and 
uniqueness in the space defined by the weighted Sobolev energy 

m pi pi 

g^oJo Jo 

when m ^ 3. In addition we show persistence of smooth solutions as long as 
the energy for to = 3 remains bounded. We do this via the method of lines, 
approximating with a discrete system of coupled pendula (a chain) for which the 
same estimates hold. 
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1. Introduction and background 

1.1. Introduction. In this paper, we explore the motion of a whip, modeled as an 
inextensible string. We prove that the partial differential equation describing this 
motion is locally well-posed in certain weighted Sobolev spaces. In addition, we are 
interested in the motion of a chain, modeled as a coupled system of n pendula, in 
the limit as n approaches infinity. We show that the motion of the chain converges 
to that of the whip. 

Although the equations of motion are well-known and have been studied by many 
authors, there are few results known about the general existence and uniqueness 
problem. Reeken |Re2j |Re3j proved local existence and uniqueness for the infinite 
string in M.^ with gravity and initial data sufficiently close (in H^^) to the vertical 
solution, but aside from this, we know of no other existence result. In the current 
paper we prove a local well-posedness theorem for arbitrary initial data for the 
finite string. 

One reason this problem is somewhat complicated is that the equation of motion 
is hyperbolic, nonlinear, nonlocal, degenerate on a spatial boundary, and possibly 
even elliptic under certain conditions. 

If 77 : M X [0, 1] M'^ describes the position ri{t, s) of the whip, then one can derive 
that the equation of motion in the absence of gravity and under the inextensibility 
constraint {rys,77s) = 1 is 

Tlu{t,s)=ds{<j{t,s)TJs{t,s)). (1.1) 

Incorporating gravity introduces some complications; to keep things as simple as 
possible, wc will neglect it. 



Equation ( 1.1 ) is a standard wave equation; however, the tension a is determined 
nonlocally, as a consequence of the inextensibility constraint, by the ordinary dif- 
ferential equation 

ass{t, s) - |?7,,(t, s)|V(i, s) = -\rjst{t, (1.2) 

With one end fixed and one end free, the boundary conditions are ri{t, 1) = and 
a{t, 0) = 0, along with the compatibility condition ds(j{t, 1) = 0. 
We use the energy 



m ^1 ^1 

e^oJo Jo 



t, s)|^ ds, 



and show that for small time we have local existence and uniqueness in the space for 
which the energy is bounded. Precisely, for any nonnegative integer m, define 



THE MOTION OF WHIPS AND CHAINS 



3 



-^m[0, 1] to be the space of functions /: [0, 1] R"^ such that 



Ik = E 



I" 

JO 



ds (1.3) 



is bounded; then E„i = \\Vt\\%^ + \\v\\n^+i- 
We prove the following result: 

Theorem 1.1. Suppose 7: [0,1] M'' and uu: [0,1] —>■ M."^ are restrictions of 
functions on [0,2] satisfying the oddness condition 7(2 — s) = —7(5) and w{2 — 
s) = —w{s), and that they have bounded weighted Sobolev norms, ||7||iV4 < 00 and 
\\w\\n^ < 00. Suppose that in addition we have 

|7'(s)P = l and (-/'is),w'is)) = for all se [0,1]. 

Then there is a T > such that there is a unique solution r/ of the equation 



(1.1) in L*([0,T],iV4[0,l]) n W^^'''i[0,T], NslO,!]) satisfying 77(0,5) = 7(5) and 



77t(0,s) = wis). 

We prove this by showing that the corresponding discrete energy 63 for the chain 
with n links is uniformly bounded for small time, independently of n. The solution 
is then a weak-* limit of the chain solutions in which converges strongly in N3 
and hence in C^. One could prove this more directly using a Galerkin method, 
but the present technique allows us to simultaneously discuss convergence of the 
discrete approximation. 

All the higher energies Em{t) can be bounded in terms of E^it), so that C"^ 
initial conditions yield solutions for short time. As a consequence, we derive 
a simple global existence criterion: if the initial conditions are C^^ functions, then 
a solution exists on [0,T] if and only if E^it) is uniformly bounded on [0,T]. 
Of course, one expects blowup of the whip equation, at least for some initial data, 
since the whole purpose of a whip is to construct the initial condition so that the 
velocity of the free end approaches infinity after a short time. See McMillen and 
Goriely |MG) for a discussion of such issues; although our model neglects some of 
the phenomena they consider, one expects that the situations are similar in many 
ways. For the heuristics of blowup in our situation, see Thess et al. |TZN| . The 
simplest blowup mechanism appears to be the closing off of a loop along the whip; 
as a loop shrinks, there appears a kink in the whip, representing blowup of both 
the curvature and the angular velocity. 

The paper is organized as follows. In Section [2| we discuss the equations for the 
whip and derive the corresponding equations for a chain with n links, in terms of 
difference operators, emphasizing the role of odd and even extensions in order to 
get the fixed endpoint conditions satisfied automatically. In Section |3] we discuss 



the solution of the tension equation ( 1.2 ) in terms of a Green function, showing that 
the tension is positive except at s = and deriving a similar result for the chain. 
We also derive sharp upper and lower estimates for the Green function. In Section 
|4]we explain why we need weighted energies, and we derive the analogues of the 
Sobolev and Poincare inequalities for weighted norms, which are used throughout 
the rest of the paper. 

In Section [5] we give estimates for the tension a in terms of 77 and rjt. For the 
norms of a we use the bounds on the Green function; for higher derivatives 
we bound the weighted Sobolev norms of a in terms of those of r] and rjt. These 
bounds are used in Section [6] to derive the main energy estimate, to bound the time 
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derivative of one energy in terms of another energy. Section [7| contains the proof of 
Theorem 1 1.1 1 Uniqueness is proved using a low-order estimate for the difference of 
two solutions. Finally in Section [8] we discuss related open problems. To make the 
paper a bit easier to read, we have moved all of the longer proofs into an Appendix. 

Victor Yudovich found several results on this problem, although he did not pub- 
lish anything on it to my knowledge. I learned of this problem from Alexander 
Shnirelman, and 1 would like to thank him for many useful discussions about it. 

1.2. Background. The study of the inextensible string is one of the oldest appli- 
cations of calculus, going back to Galileo, and yet it is still being studied to this day. 
One is especially concerned about kinks in the solution and what the appropriate 
jump conditions should be; authors such as O'Reilly and Varadi [OVj . Serre [Se] . 
and Recken [Relj have discussed these issues in detail from differing points of view. 

The first problem to be studied was finding the shape of a hanging chain, first 
solved incorrectly by Galileo and then correctly by Leibniz and Bernoulli, one of the 
first major applications of the calculus of variations. The shape of small-magnitude 
vibrations of a chain hanging straight down (in a linear approximation) goes back 
to the Bernoullis and Euler [^, and is taught in textbooks today as an example 
of Bessel functions; see Johnson [J] and Schagerl-Berger [SB] for related problems. 
Kolodner jKo| . Dickey |D1| . Luning- Perry |LPj . and Allen-Schmidt |ASj studied the 
problem of a uniformly rotating inextensible string, one of the few other problems 
that can be solved more or less exactly. 

Burchard and Thomas [BTj obtained a local well-posedness result for the related 
problem of inextensible elastica, in which there is a potential energy term reflecting 
a resistance to bending; however it is not clear whether the solutions are preserved 
in the limit as the potential term goes to zero, so this result does not help in the 
present situation. 

Many authors have studied the problem of a vertically folded chain falling from 
rest; this is a classical problem that appears in several textbooks ( jAnj . jD2j . [H] . 
and |Ros| ) . In recent years the problem has been debated in the physics literature, in 
particular the issue of whether energy is conserved and whether the tip of the chain 
falls at an acceleration equal to gravity or faster ( |Calj [CalMar] |CapMaz| |dSR| 
[HHR] [ffl] [OV] [SSST] [TP] [TPG] [ST]). See Wong-Yasui [WY] or McMiUen [M] 
for a good survey of the literature. 

McMillen and Goriely ( [GM] and [MS]) studied a tapered whip theoretically, 
numerically, and experimentally, showing that the crack comes not from the tip but 
rather from a loop that straightens itself out. They use a different model, however, 
in which the tension depends locally on the configuration. Thcss et al. |TZN] 
studied the blowup problem for the closed inextensible string, especially as a model 
of the blowup problem for the Euler equations for a 3D ideal fluid. They found 
evidence of blowup from loops closing off, showing numerically that sup^|?7st| ~ j^^^ 
and supg|77ss| ~ (T-f)^/^ ' '^^^^^re T is the blowup time. 

2. The basic equations 

In this section, we present the equations for both whips and chains, assuming 
no external forces. Our boundary conditions come from the assumption that one 
end of the whip or chain is held fixed at the origin, while the other end is free. We 
describe the whip as a function r] : [0, T] x [0, 1] ^ E*^, and describe the chain as a 
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sequence of functions rjk'- [0, T] ^ M'' for 1 ^ fc ^ n + 1. Our formulas simplify if 
we assume the fixed point occurs at s = 1, i.e., ry(i, 1) = for all t; for the chain, 
we assume ?7„_|_i(i) = for all t. 

2.1. The whip equations. We will just present the equations here with a sketch 
of the derivation; the reader may refer to ^ for a detailed derivation and discussion. 
Schagerl et al. jSSSTj and Thess et al. [TZNj also present derivations from minimum 
principles: the basic idea is to find a critical point of the action ^^\7]t{t, s)P ds dt 
subject to the constraint \r]s(t,s)\'^ = 1. 

A variation ( must satisfy ({t, 1) = and (r]s{t, s),dsC{t, s)) = 0, and if 77 is a 
critical point of the action, then '^^(rjt, dtQ ds dt = for all such C- Integrating 
by parts, we conclude that a critical 77 must satisfy the evolution equation 

Tju{t,s) = ds{<j{t,s)rjs{t,s)) v{t,l) = 0, (2.1) 

for some function cr, where a{t, 0) = 0. 

Differentiating |?7sp = 1 twice with respect to t, we find that a is determined 
by the following boundary-value problem for an ordinary differential equation (for 
each fixed t): 

assit,s)-\r,ss{t,s)\''ait,s) = -\r,stit,s)\^, a(i,0) = 0, (T,(i,l) = 0. (2.2) 

The boundary conditions are compatible with the evolution equation as long as t] 
can be extended to an odd function through s = 1; in that case a can be extended 
to an even function through s = 1, which is where we get the extra boundary 
condition (Ts(i, 1) = 0. See Figure [l] Oddness and evenness give us the correct 
boundary conditions for all higher derivatives of rj and cr at s = 1, which is crucial 
for the a priori estimates. Furthermore there is a discrete analogue of oddness 
and evenness for the chain which both simplifies the equations and helps greatly in 
defining the higher discrete energies. 







s = 




s=l . 




'■••...5 = 2 





Figure 1. The free end of the curve is at s = 0, while the fixed 
end is at s = 1. We imagine the curve extending smoothly through 
the origin to s = 2 through an odd reflection such that r]{s) = 
—77(1 — s). Under such an extension, the tension extends to a 
smooth function satisfying (t(s) = <t{1 — s). 
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A sometimes simpler way of dealing with the constraint |?7s| = 1 is to consider 
i]g as a curve on the unit sphere in S"^^ ^ . For simplicity we will assume d = 2 when 
doing this, although the technique works in spherical coordinates in any dimension. 
We write 

T]sit,s) = ( cos 6'(t,s), sin 6i(t,s)); (2.3) 
a straightforward computation verifies that ( |2.1[ ) becomes 

dtt{t, s) = a{t, s)9,s{t, s) + 2as{t, s)9sit, s), (2.4) 



while (2.2 1 becomes 

ass{t, s) - Osit, sfa{t, s) = -Otit, sf. (2.5) 
The fact that 77 is odd through s = 1 forces us to have 9 even through s = 1, 



so the boundary condition on (2.4) is 9s{t,l) = 0. We could work out all the 



estimates directly in terms of the system (2.4)-(2.5), but the discrete versions of 



these equations are substantially more complicated than the discrete versions of 



(12JJ)-(12^, even when d = 2. 

If a-{t,s) is strictly positive for < s ^ 1, then equation ( |2.1[ ) is a hyperbolic 
equation with a parabolic degeneracy at s = (since we must have a{t, 0) = 0). As 
such, the only condition necessary to impose at s = is that r]{t, 0) remain finite. 



We point out that equation (2.1 1 cannot be an ordinary differential equation on 
any infinite-dimensional Sobolev manifold: the right side is obviously an unbounded 
operator even in the simplest case. Hence we cannot hope to prove existence and 
uniqueness using the techniques of Picard iteration on an infinite-dimensional space, 
as in Ebin-Marsden jEMj . Instead we will work directly with the partial differential 
equation using energy estimates. 

2.2. The chain equations. We now derive the equations for the finite model, 
consisting of {n + 1) particles in R'^, each of mass ^, one of which is held fixed. 
The particles are assumed to be joined by rigid links of length ^, whose mass is 
negligible. The position of the fc*^ particle is ?7fc(i) for 1 ^ fc ^ n -I- 1; we assume 
the fixed end is the {n + 1)*'* particle, so that r]n+i{t) = for all time[^ The 
configuration space is thus homeomorphic to {S'^^^)'^, and is naturally embedded 



m A.'^". 



The kinetic energy in M. is 



2n ^ 



1 " 

;;;i:\v>^\'- (2-6) 



In addition the constraints are given by 

hk{r]i:---,Vn)=^\Vk+i-mf = l^k^n. (2.7) 

Stationary points of the constrained action satisfy the equations of motion ry^, = 
'^^"^j ^nk^j ^'^^ some Lagrange multipliers aj. More explicitly, we have 

rjfc = n^akiVk+i - Vk) - TT?Uk-i{m - Vk-i) (2.8) 

for 1 < fc ^ n. The scaling by is chosen so that ak{t) converges to a function 
<j{t,s) as n —>■ CO. The numbers a physically represent the tensions in each link. 
We set (To = so the same equation is valid when k = 1. 

^It might seem more natural to assume »?o(i) ~ 0, but our choice makes the tensions (7^ 
proportional to k/n rather than (n — k)/n, simplifying some formulas. 
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The constraint equations determine the a. Differentiating (2.7) twice with re- 
spect to time and using (|2.8[), we get 



+ n'^cTk-iiVk - Vk-i,Vk+i - Vk) (2.9) 



for 1 ^ fc < n (again using (Tq = 0), while for k = n we get (using rjn+i = 0) 



We note that if 



1 " 

?7fe(i) = - - 2 n) '^k{t) = (j{t, 



(2.10) 



(2.11) 



where 77; E x [0, 1] 
formal limit of (2. 



and (T : 



is pT)) and the formal limit of ([2^1 is 
1 for all k 



discretization of 77 ensures that KV+r/)^ 
j. We will refine this in Section [7.1| ) If ri{t,l) 
choice also gives 77„+i(i) = and cro{t) = 0, 



[0, 1] ^ M are C°^', then as n — >■ 00, the 

(Note that this 
since \ris{t, -)| = 1 for all 
and a{t,0) = 0, then this 
as desired. Hence the chain equations 



(2.8) and (2.9 1 form a discretization of the whip equations (2.1) and (2.2) which 



conserves energy as well as preserving the geometry. 

The analysis of the chain equations becomes much simpler if we can avoid using 
separate equations for the boundary terms. An easy way to do this is to extend rjk 
and cTfe beyond fc = n by demanding that ry be odd through k = n + 1 and that a 
be even, which is exactly what we had to do for the whip in Section |2.1[ So for 
fc ^ n + 1 we set 



11k = — ?72Ti+2-fc, 
Cfc = (T2Ti+l-fc- 



(2.12) 



Then it is easy to see that the evolution equation (2.8) still holds for the fixed point 
at fc = n + 1 and that (2.9) for k = n yields the tension boundary condition ( |2.10 ). 

A further simplification comes from using difference operators. (See for example 
[LLj .) First recall that for a sequence / defined on some subset of Z, the (forward) 
shift operator E is given by {Ef)k = fk+i- The backward shift is denoted by 
E^^, so that {E^^f)k = /fc-i, and powers of E signify composition. We define the 
(forward) difference operator V+ by 



(V+/)/c = n[fk+i - /fc], 



(2.13) 



so that if / denotes the identity operator, then V+ = n{E — I). It is also sometimes 
convenient to work with the backward difference operator V_, defined by (V_/)fc = 
n[fk - fk -i], so that V- = E^'^\/+ = n{I — E^^). In this notatiorj^ equations ( |2.8[ ) 
and (|2.9|) become 



i) = V_(crV+ry), 
<V+r;,V_V+(aV+r;)> = -|V+77p, 



(2.14) 
(2.15) 



The more usual finite-difference notation is A for the forward difference and V for the back- 
ward difference; we use V-i- and V_ instead to avoid confusion with the Laplacian on smooth 
functions, and since our rescaled version is not standard. We prefer the rescaling since if the 
sequence /j, converges to a smooth function f(s) as n — > c», then (V-i-/)^ converges to f'{s). 
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where the equations are vahd when any subscript 1 ^ k ^ n is placed on aU the 
terms simuhaneously. We can thus write all the discrete equations without specific 
reference to subscripts, which simplifies the notation. 

The following formulas will be useful when working with difference operators and 
sums: both follow from the simplest product formula V+(/.g) = 5V+/ + EfW+g. 



- y! 9k^+fk = - -Yj fk^-9k + fn9n 



/o5o 



(2.16) 



(2.17) 



fc=0 



We can rewrite (2.15) in a more useful form, solving for the second difference 
V_V+cr in terms of everything else, using |V+ryp = 1 to simplify the terms. We 
obtain 

V_V+a = ^iV^+ryp + :^|V_V+77p - (2.18) 



and the resemblance to the continuous version (2.2 1 is obvious. 



3. The Green function for the tension 



At each fixed time, equation (2.2 1 is a linear nonhomogeneous ordinary differen- 
tial equation for the tension a . Hence there is a Green function G'(t, s, x) depending 
on such that 



f G( 
Jo 



a{t,s)= G(t,s,x)\r\tx(t,x)\ dx. 



Similarly, equation (2.18) can be thought of as a linear nonhomogeneous matrix 



equation for a, for which the solution takes the analogous form 

1 " 

for some "discrete Green function" Gkj- Naturally we expect that if ^ ^ s and 
—>■ X, then Gk„j„it) —>■ G(t, s, x) as n —>■ 00; this can be proved as a consequence 
of our general convergence result for rj. Our goal in this section is to establish 
properties of these Green functions. In particular we establish that the Green 
function is always nonnegative for a whip in a sufficiently smooth configuration, 
while the Green function is nonnegative for a chain as long as all the angles between 
links are obtuse. Furthermore we want to establish upper and lower bounds for the 
ratios and '^''jj*^ , where Sk = in order to be able to compare the norms 

weighted by powers of a{t, s) to the norms weighted by powers of s. 

3.1. Basic properties of the Green functions. First we discuss the solution op- 
erator of the whip tension. To keep the notation relatively simple, we will suppress 
the time dependence. 



Proposition 3.1. For any fixed time t, the solution a{s) of (2.2) is given by 

a{s) = f G{s,x)\T]tAx)\^dx, (3.1) 
Jo 
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where G is the Green function given by 

Gss{s,x)-\'nss{s)\''G{s,x) = -5{s-x), G(0,x) = 0, G,(l,a;) = 0. (3.2) 

The Green function is symmetric, i.e., G{s,x) = G{x,s). It satisfies G{s,x) > 
whenever x > and < s ^ 1. In addition if < x < 1, we have Gs(s, a;) > for 
< s < X and Gs(s, x) ^ for x < s < 1. 

Proof. The existence of the Green function and the symmetry property G{s,x) = 
G{x, s) is a well-known result of the general theory for second-order equations with 
homogeneous boundary conditions. See for example Courant-Hilbert |CH] . 

To prove the other statements, we first show that G{x, x) > for any x 6 (0, 1). 
For any fixed x, multiplying (3.2 1 by G{s,x), integrating from s = to s = 1, and 
using integration by parts with the homogeneous boundary conditions shows that 



G{x, x) 



f Gs{s,x)^ds+ r 

Jo Jo 



r]ssis)\'^G{s,x)^ ds, 



which forces G{x, x) ^ 0. Because of the jump condition 

lim Gs{s,x) — lim Gs(s,a;) = l, 

s— >2;~ s— >a:+ 

we cannot have Gs{s, x) identically zero if < x < 1. So G(x, x) > if < a: < 1. 
It is then easy to prove the other statements in the intervals (0, x) and (x, 1) using 
the boundary conditions. □ 

Now let us do the same for the tension operator for the chain. The equation 
(2.9), or the more elegant version (2.18), makes clear that the vector (cri, . . . , (t„) of 



tensions comes from inverting a tridiagonal matrix. Since this is one of the easiest 
matrix types to invert, we get a relatively explicit formula for the solution, which 
will be useful in constructing estimates on the maximum and minimum tension. 



Proposition 3.2. The solution of the constraint equations (2.15) is 



-Vol 



(3.3) 



where the discrete Green function G^j is constructed by 



G 



min{j,fc} 
1 PijPik 



whe 



J-1 

n 



a, = <V+7/,+i,V+77,> (3.4) 



and P satisfies the recursion 

Pn = 1, ^i^ 



for 1 ^ i ^ n — 1. 



In (3.4 1 we use the convention that the empty product when j = i is 1. 



(3.5) 



The tensions ak are positive for every nontrivial choice of V+rj if and only if 
ai > for every i. 



Proof. The system (2.9) and (2.101 is of the form Aa = w, where A is a symmetric 



nonnegative diagonally-dominant tridiagonal matrix and w is the vector of angular 
velocities Wi = |V+?7ip. There are several standard algorithms for inverting such a 



matrix; the formula (3.3 1 is given in the review paper of Meurant [M] . 

Clearly af ^ 1 for all i, so that inductively we have 1 ^ /3i ^ 2 for all i. Hence if 
all Ui are positive, then all Pij are positive and hence all terms Gkj are positive for 



10 STEPHEN C. PRESTON 

1 ^ j, k ^ n. Thus if any V+T^j is nonzercj^ then (3.3) says that all ak are positive 
for 1 ^ A: ^ n. It is easy to see that if a, for some i, then there is some choice 
of V-)-?7 so that some cr is nonpositive. □ 




Figure 2. In the top row we plot chains with constant angles 
between consecutive segments, with n = 4, n = 8, and n = 56 
links. In the bottom row we plot the corresponding discrete Green 
function Gkj as a function of ^, evaluated at j = for each 
configuration. Notice that when n = A the angles in the chain 
are acute, which is what allows the tension to become negative in 
that case. Also notice that as n ^ oo, the discrete Green function 
approaches the Green function for the differential equation. 



3.2. Upper and lovifer bounds for the Green functions. Proposition [XT] im- 
plies that if < a; ^ 1, then G{s,x)/s is a positive function of s on [0, 1], since 
lims->oG{s,x)/s = Gs{0,x) > 0. We now want to know exactly how large or 
small this positive function can be; ultimately our interest will be in the quanti- 
ties supjgj-Q ^ cr(s)/s and iufsgjo^i] a{s)/s, which are completely determined by the 
bounds on the Green function. We are especially interested in the discrete ana- 
logues, maxi^fc^n nGkj/k and mini^fc^„ nGkj/k. We end up with the same upper 
bound in both cases, which is relatively easy to prove, while the lower bound is 
much more complicated and necessarily weaker in the discrete case. 
First we establish the upper bound. 

Proposition 3.3. If rj^g is smooth, then the Green function G{s,x) defined by 
Proposition \3. 1\ satisfies the following bounds. 

G(s x) 

sup \Gsis,x)\ ^ 1, and sup ' — ^ 1. (3.6) 

O^s.x^l 0^s,a;^l S 

Furthermore, suppose Gkj, rjk, ctk, o.nd (3k are as defined in Proposition \3.S\ and 
that ak = {V_|_77fc_|_i, V_|_?7fe) ^ for all k, so that Gkj ^ for all j,k. Then if 



'Of course, the only way every V+t^j is zero is if the chain is stationary, since »7n+i = always. 
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(V_,iG')fcj denotes the partial difference (V 
convention Gqj = 0, then 

nGkj 



iG)kj = n{Gkj - Gk-i,j), using the 



|(V_,iG')fcj| ^ 1 and 



k 



^ 1 for all 1 ^ j,k ^ n. 



(3.7) 



Proof. The proof of (3.6) is easy: by Proposition 3.1 the partial derivative Gs(s, x) 
is positive for s < x and nonpositive for s > x, and jumps by — 1 at s = x. 
Since Gss ^ whenever s ¥= x, we know Gg is increasing on each interval. We 
therefore must have < Gs{s,x) ^ 1 for s < x and —1 < Gs{s,x) ^ for 
a; < s ^ 1; either way, |Gs(s,x)| ^ 1. Then using the fact that G{0,x) = 0, we 
have ^ G{s,x) = drG{r, x) dr ^ dr = s, which yields (3.6). 

The proof of 
from Proposition 



(3.7) is more complicated, but uses the same basic ideas. First, 
we know that Gkj ^ for all k and j since every ^ 0. 
Assume first that j ^ n. Then rewriting (2.9)-(2.10), we see that the discrete 
Green function satisfies the equation 

1 



ctk+iGk+ij 
1. while 



2Gkj + ak-iGk- 



-Sk 



k] 



for 1 ^ fc < n 

— Gnj + a,i_iG„_i J- = 

Since ak = (V+Ty^, V+ryfc+i) with |V+77fc| = 1, we have ak ^ 1. Thus for fc # j we 
can easily see the second partial difference satisfies (V+_i V_,iG)fcj ^ 0. Since the 
second partial differences are nonnegative except at the diagonal, the first partial 
differences are increasing except at the diagonal, i.e.. 



When k 



(V_,iG)fc+ij - (V-,iG)fc, ^ for aU k ^ j. 
■ j we can check that 

(V_,iG),+i,, - (V_,iG)„- 



(3.8) 



(3.9) 



Now look at the endpoint terms: at the left endpoint, we know (V_,iG)ij 



Gij ^ 0. At the right endpoint, if j ^ rt then we have —Gnj + Gn 



nL _ 

an-i)Gn-i,j ^ 0, so that (V_aG),y 
conclude that if j n then 

<; 



(1 



^ 0. Thus combining (3.8) and (3.9), we 



sc (V_,iG)ij • • • sc (V_,iG)j, ^ 1 + (V_aG),+ij sS • • • 1 + (V_aG)„, 

Hence we must have |(V__iG)fcj| ^ 1 for all fc, as long as j n. 
If j = n, the situation is slightly different; in that case we get 

(V_,iG)i„ (V_,iG)2„ ^ • • • (V_,iG)„_i^„ s; (V_^iG), 

so that |(V_^iG)fej| ^ 1 even if j = n. This completes the proof of (3.7|. 



s; 1. 



s; 1, 



□ 



Remark 3.4. Unfortunately we cannot bound '^^^^^^ from above. If we denote by 
Go(s, x) the Green function when \riss\ = 0, then we easily compute that Go(s, x) = 
min{s, x}, so that '"^"^^'^^ = min{i,i} is unbounded on [0,1] x [0,1]. Note that 
by the Sturm comparison theorem, we have that G{s,x) ^ Gq{s,x) for any Green 
function satisfying (3.2). However it is easy to see that for any ^ p ^ 1, we have 

G{s,x) ^ Gq{s,x) 



^ 1. This will be useful in the proof of Theorem 



7.7 



sPx'--P sPx^ 

It is easy to check that the discrete Green function satisfies the same inequality, 
\Gkj\ ^ ^mm{j,k}, using formula (3.4) and the fact that \ai\ ^ 1 and A ^ 1 for 
all i. In fact this bound is valid even if not all Gkj are positive. 
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Now we establish the lower bound. This is the only time in the paper where we 
get a weaker result for the chain than for the whip; the reason is that we need to 
make strong assumptions in order to prevent sharp kinks in the chain, to ensure 
nonnegative tension. Smoothness of the whip, on the other hand, ensures that the 
tension in the whip is nonnegative automatically. 



Proposition 3.5. Suppose Gkj, Tjk, at, and fik are defined as in Proposition 3.2 
Assume the rjk are such that, for some v 6 (0, ^^], we have 

^3/2 

-Tj7^\ylVk\^ ^ V for alll k (3.10) 
Then for all 1 ^ j, fc ^ n we have 

^ e-^". (3.11) 



jk 



If G solves (3.2) and rjss is a smooth function, then we have 



Proof. The two estimates are proved in slightly different ways, but the main point 
for both estimates is to show that the minimum is attained at the off-diagonal 



corner, then estimate this value either using the direct formula (3.4 1 (for the chain) 



or through a substitution (for the whip). The full proof is in Appendix A.l □ 



The assumption ( 3.10 ) for the discrete case is much stronger than the assumption 



JO 

to ensure every > in order to get all tensions positive (by Proposition 3.2 1, 



s\r]ss\'^ ds < 00 for the continuous case, but such a pointwise bound is necessary 



even when k = 1. The exponent | is important: the exponent 1 would work to 
prove the estimate, but we cannot prove that such an estimate actually holds for 
all values of t; the exponent 2 is not enough to get a lower bound for the tension. 



Remark 3.6. Note that we could easily get a stronger estimate than (3.12) if we 
simply assumed an upper bound on \r]ss\, using the Sturm-Liouville comparison 
theorem. However, we prefer the weaker assumption that sjTyssp ds < oo, since it 
allows for the possibility of the curvature at the free end of the whip approaching 
infinity (a possibility not precluded by the equations due to the degeneracy there). 
Even if the weighted energy is finite — the condition under which we will prove 
local existence — we will not necessarily have \r]ss\ bounded on [0, 1]; an example is 
when 6{s) = s'' for some q€ (^,1), using the spherical representation (2.3). See 
Example |5.3| for details. 

4. Weighted Sobolev norms 

4.1. Motivation. In order to demonstrate existence and uniqueness, we want to 
apply the usual technique of energy estimates in Sobolev spaces. By showing that 



we have sequences of solutions of the chain equations (2.14 )-( 2.15 1 for which the 
energy is uniformly bounded, we can extract a convergent subsequence to establish 
existence; the same sort of energy estimates can also be used to establish uniqueness. 
Several issues arise to complicate this strategy. 
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Ordinarily for a wave equation like (2.1 ), one would try to bound an energy like 

Eo = f \vt\^ + <j\Vs\'' ds 



by computing its derivative and using Gronwall's lemma. It is easy to compute 
that 



dEp 
dt 







at{t,s) 



^ ds ^ [ sup 



Eoit), 



using the boundary conditions cr(i, 0) = and rit{t, 1) = 0. Unfortunately we cannot 
bound at or even a in terms only of Eq. Indeed, it is hard to even make sense of 
equation (2.2) unless both rjst and ij^s are in L^, which means we have to consider 
higher energies. 

Here a complication arises. The usual approach would be to consider an energy 
like 

Fi = Eo+ r \rist\^ + alijssl"^ ds. 
Jo 



Its derivative is, using (2.1) 
dPi 



dt 



< 



(Jt{t,s) 



Flit) + 2 \ <Js(v. 
Jo 



t,riss)ds. 



Here the boundary term vanishes since (T{t,0) = and rjssit,!) = (recall we 
assume 77 extends to an odd function through s = 1). Furthermore since \r]s\'^ = 1, 
we have (rjs, rjst) = 0. The problem is that if we want to get the right side in terms 
of -Fi alone, we need to use the Cauchy-Schwarz inequality to get 



J" 

Jo 



as{r]st,ilss)ds 



VosSssSl ^<7[t,s) 



hut the right side is not bounded. We always have cr(i,0) = 0, while we will 
generally not have crs(i,0) = 0. 

Instead we want an energy for which the integration by parts cancels out this 
highest-order remainder. The only such quantity of the form v4|77stp + -B|ctssP ds 
for which this works is 



Ei = Eo+ { <j\rist\^ + a'^\riss\'^ ds. 
Jo 



With such a choice we get 



dEi 
"df 



2 sup 



O^ssSl 



crtit,s) 
(j{t,s) 



Elit), 



which we can manage once we understand how a and at behave. The same phe- 
nomenon continues for the higher energies as well, which motivates us to define 



m pi 

Em = T. a'\8%\' + a'^'\di^'v\'d.s. 
e^o-'o 



(4.1) 
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With this definition, we have 



dt 



L Jo Jo 



Y,2CY) lla'di+^-^a(d%,dl+^r,)ds 



i=0 



(4.2) 



with the three remaining terms integrating to give cr^+^<5f77t, 5f+^77)|g=J = due 
to oddness of r] through s = 1. 

Our primary goal will be to bound (4.2 1 in terms of the energies Em- More 
specifically, using the fact that a{t, s) degenerates like s near s = 0, we want to get 
bounds in terms of the simpler weighted energies Em defined by 



En 



ds. 



(4.3) 



To do this, we will need several estimates. So our first goal is establishing basic 
Sobolev-type and Wirtinger-type inequalities for such weighted norms. In addition 
we need to show the energies (4.3) are equivalent to the tension-dependent energies 
(4.1), which means we have to bound sup^ cr(i, s)/s and iidg a{t, s)/s away from 
zero. (The constants in these bounds will also turn out to depend on the energies 
(4.3).) Most of the work for this was done in Section |3] 

4.2. Definitions and properties of weighted seminorms. First let us define 
the weighted Sobolev and supremum seminorms we need. 

Definition 4.1. The weighted Sobolev seminorm of a function /: [0, 1] M'' is 

defined by 



Iff 

\J Mr. 71 



Jo 



/(")(s)f ds. 



The weighted supremum seminorm of / is 



2 

r,m 



sup snp"'>{s)\ 



We want to define a discrete analogue of each of these, for a sequence {/i, 
with values in M''. For this purpose, it is convenient to set 



T{k + r) 



for k e {1, . . . , n} and for any real r > — 1, 



(4.4) 
(4.5) 

) /n} 

(4.6) 



where F is the usual gamma function satisfying F(a; + 1) = xT{x) for x > and 
F(fc) = (fc — 1)! for k a natural number. These are rising factorials, which are more 
convenient for our purposes than the falling factorials typically used in difference 
equations; cither is much more convenient in studying difference equations than 
simply using the powers {^Y; see |LL| . Clearly if fc„ is a sequence such that 



lim„_ 



s, then we have lim„_ 



,(0 



Crete analogues of (4.4|-(4.5) by 



Recalling the definition (2.13) of the difference operator V+, we define our dis- 

-j n—m 

- E 4^^|VV"A-P (4.7) 



I/I 



2 

r.m 
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and 



2 



max 4^)|V™A.P. 

l^k^n—m 



(4.8) 



We use the same notation for both norms to emphasize the analogy; every estimate 
we prove for the discrete seminorms (|4.7|)-(|4.8|) wiU have constants independent of 



n, so that we get corresponding estimates for the smooth seminorms (4.4|-(4.5). 
Clearly if / is smooth and we define fk = /(^) for each n, then lim„^oc' |1 /« II r.m = 
||/||r,m and lim„^oc, |||/n||| = |||/|||. 

Now let us describe the main estimates. For two norms ||-||i and ||-||2 on func- 
tions, we use the notation ||/||i < ||/||2 to mean ||/|| ^ C||/|| for some constant C 
independent of /. If / is instead a sequence, then this notation will imply that C 
is also independent of n. 

For unweighted norms of smooth functions, we have the Wirtinger inequality 



ri/(^ 

Jo 



■ds < 



We also have the Sobolev inequality 



sup |/(s) 



< 



f f{s)ds \ r \f{s)fds 

Jo JO 

C\f{s)\^ds+ C\f'{s)\^ds. 

JO JO 



(4.9) 



Our weighted versions of each are as follows. 

Theorem 4.2. Let f : [0, 1] R'^ be C*. Then for any r > the norms 
and (4.5) satisfy the weighted inequalities 

ll/llr-l,m ^ ll/llr,m + ll/llr+l,m+l 
/ 1 1 r,7n 11/11?' 



2 

r,rn 



+ and (4.10) 

+ ii/ii;+i,m+i- (4.11) 

If in addition we have /(™)(1) = 0, then these inequalities can be simplified to 

\\f\\l-i,n ^ ll/ll'+i,™+i ^nd (4.12) 
Ill/Ill'™ ^ ll/ll'+i,™+i- (4.13) 
// / is instead a sequence {/i,--- ,/«} with values in M."^, then the inequalities 



(4.10 )-(4.11 ) also hold if the norms are interpreted as (4.7) and (4.8|, while the 

0. 



inequalities (4.12|-(4.13) hold if f^- 



Am) 



Proof. The continuous version of this inequality appears in Adams- Fournier jAFj . 
We prove the discrete version in Appendix |A.2[ from which the continuous version 
follows in the limit. □ 

Remark 4.3. The example f{s) = arcsinh (Ins) demonstrates that the inequalities 
(4.12) and (4.13) cannot be extended to r = 0: in that case we have /(I) = 0, 
CO, and s\f'{s)\'^ ds < co, while l\f{s)\'^ ds and sup^g^ |/(a:) p 
are both infinite. In particular there cannot be constants for the discrete versions 
that are independent of n when r = 0. 



The important thing about (4.11 1 is that by Remark 4.3 the estimate only works 
when r > 0. Hence in any computation where a supremum norm is required, we 
will want a positive power of s attached to be able to use this result. This will 
show up when we need to estimate weighted Sobolev norms of products of three 
functions: we want to pull out a supremum of one and use Cauchy-Schwarz on the 
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rest, and we will need a little extra weighting in some cases. Of course, we could 
use the usual Sobolev inequality (4.9 1 to get 

|||/|||0,0<||/||0,0+||/||04^ll/l|0,0+|l/|ll,l + 11/112,2, 

but requiring two extra derivatives rather than one is usually not worthwhile (except 
once in the proof of Theorem 7.7 1. 

Frequently our weighting in discrete norms will be slightly off (for example, we 



may want to replace sl,^''''^Vsl'^ with s^^^ 



k ' 



or we may want to replace sS?'^ with 



for some j). In the continuous case these formulas are trivial, but in the 
discrete case, bounds such as these come from properties of the gamma function 
(in particular the fact that the gamma function is log-convex by the Bohr-MoUerup 
theorem). The constants will never be important; what will matter is that they are 
independent of k and n. The following estimates are easy to prove. 

Proposition 4.4. Let n e N and let k e {1, . . . , n}. Let p and q be positive real 
numbers. 

Then the weight function si^'^ = ^^p"tf] satisfies the following inequalities: 



k nPr{k) 

ip+q) 



rip + q + i) 



„(p) 



We also have ^ s^^^. ^ r(j+i)r(p+i) 



Jp) 



Tip+mq+1)^'' • 
^+i)r(" +1) Sfc'^ /o*" ^.TT-y nonnegative integer j . 



(4.14) 



Proposition |4.4| also gives the following corollary, which is the most useful tool 
we have for estimating norms of products. To get the higher-difference norms of 
products, we will use the product rule (2.16) for differences together with these 
formulas. The proof is trivial. 

Corollary 4.5. Suppose (/i, . . . , /„) and (gi, . . . ,gn) are sequences of real numbers, 
and let p and q be nonnegative real numbers. Then 

WfgWM ^ lll/lilp,oll5ll^,o (4.15) 



P+9,0 



< 



2 llloIlP 
p,Ollli/lllrj,0- 



(4.16) 



and the other in 



The formulas are also valid if one of the sequences is in 
if both are in M'' and we use {fjg} or \f\\g\ instead. 

Estimates (4.15) and (4.16) are also valid if f and g are smooth functions with 
the norms interpreted as (4.4) and (|4.5|). 



Remark 4.6. Typically we will extend the seminorms (4.7) and (4.8 1 when used 
for r] and a; since we have ao = and ?7n+i = 0, it is more convenient to modify 
the definitions to 

-J n— m+l 



n— m 



fe=0 



This does not affect any of the estimates, but it allows us to incorporate the endpoint 
information. This is convenient for example to interpret (3.3) in terms of WfjWl i in 
Lemma 15.11 
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4.3. Weighted energy norms. We have already defined the weighted energy (4.1 ) 
and (4.31 for a whip. Now we want to define corresponding discrete energy for the 
chain. The definitions are made much easier if we use the odd extension (2.12) of 
f] to define the differences V™?7fc beyond k 
with (4.6), we define 



(r) 



n — m + 1. Furthermore, by analogy 

(4.17) 



fc+r-l 

Y\ f j for any integer r ^ 0. 

j^k 

Our time-independent energy will be 

, m ri-[i/2\ 



(4.18) 



(4.19) 



=0 fc=l 

while the time-dependent energy is 

^ m ri-[e/2\ 
n 

£=0 fe=l 

Recall that we need to use the time-dependent cr-weighted quantities to compute 
the time-derivative of energy in order to get some cancellation, while only time- 
independent energies are useful for constructing topologies and relating distinct 
norms. 

Clearly if we have sequences r\n and (t„ defined for each n e N as in (2.11 ), then 
Em\r\\ = lim e,n[iln] and E,n[ri] = fim e.m[??n]- 

n— n— 

So any estimate we obtain on the chain energies Cm and will become an a priori 
estimate on the corresponding whip energies and Em- 
Note that we have 

m m + 1 

en.^J]MlE+ J]Mli, (4.20) 

in terms of the discrete weighted seminorms ( |4.7[ )p| 

It is also convenient to observe that the lowest-level energy eo is constant in time. 
Since it will be useful later in Lemma 5.2 we separate the terms and define 



Wo 



1 " 
7. Xi 



ml,. 



1 " 

n 2 



Sfc|V+77fe| 



(4.21) 



fe=i 



Lemma 4.7. Ifrj satisfies (2.14) with |V+7y| = 1, then both uq andvo are constant 
in time. 

Proof. The fact that |V+r/| = 1 implies that = ^ + As a corollary, every 
energy given by (4.18) satisfies ^ 



For uq, we just compute 
duQ 2 



Yi('nk,m) = - 2<?7fc, v_((TV+r;)fc>. 



k = l 



We would have equality if the sums over k went from k = I to k = n — £ rather than 
k = n — [^/2J. The reason the sums in ( |4.18[ l and l |4.19| contain a few extra terms in the sums is 
in order to make the derivative estimate of Theorem |6.1| simplcr: with this definition the endpoint 
terms of the discrete energy derivative always vanish, as they did in ||4.2|. 
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Using the summation by parts formula (2.17) along with (V-t;, V-?)) = and the 
endpoint conditions cfq = and tjn+i = 0: it is easy to show this sum vanishes. □ 



For a smooth solution of ( |2.1[ )-( |2.2[ ), we clearly have that the analogous quanti- 
ties Uq = ^^\r]t{t, s)p ds and Vo = s)P ds satisfy Vq = ^ and Uq is constant 
in time. Thus Eq is also constant in time. 

Our primary use of Theorem |4.2| will be the following formulas, which follow 



easily from (4.20) 



Lemma 4.8. For any i ^ and ^ j < i, we have 
\\v\\t-j,i^e,+j-i, 



\\v\\\i-j,i ~ ^i+ji 
lll'7llll/2,i 



< 



e2i+l- 



(4.22) 
(4.23) 
(4.24) 



Remark 4.9. We need the extra power of 1/2 in estimate (4.24), since (4.11 
not valid when r = 0. This is important once at the end of the proof of Lemma 



IS 



5.4 



and once at the end of the proof of Theorem 6.1 



5. Bounds for the tension in terms of the energy 
Before bounding the energy itself, we first want bounds for the tension a given 



by either (2.2) or (2.15). To compare the energies Em and Em, we want upper and 



lower bounds for a/s. In addition, to compute the time derivative of Em, we need 
to know a bound for at/s, by formula (4.2). 

For a smooth solution (?7,cr) of (2.1)-(2.2), we define quantities A, B, and C by 
the formulas 

Ait)= sup \asit,s)\, B{t)= sup C{t) = sup |cr^f(i,s)|. (5.1) 

Observe that since (j(t,Q) = 0, we have 



\ait,s)\ 



Jo 



dx 



^ s sup \ax{t,x) 



so that we have 



sup ^ ' sc Ait), and similarly sup ' ^ Git). 



(5.2) 



Generally the bounds (|5.2[) will be much more useful to us, although occasionally 



we will need the actual definition (5.1 ) 



Similarly, we define the discrete analogues of the quantities (5.1); as with the 



energy, we use upper-case and lower-case for norms of the whip or chain respectively. 
Recall that = ^, while our convention is that cro(i) = 0, and recall the definition 
(V_cr)fc = 7^(c^A; — CTfc-i). We therefore set 



ait) = max |(V_cr)fc(t)|, b{t) 



max 

l5=fcs=n crfc(t) 



c{t) = max \{\/-a)kit)\. (5.3) 



As above, the fact that cro(0 ~ means we can write 

k 



Wk{t)\ 



1 fc / 

-I](V-fT),(t) ^ - max|(V_a 

.7 = 1 ^ 
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to obtain 



CTfc 

max — ^ a 



and 



max ^ c, 

l^fc^n Sk 



(5.4) 



Lemma 5.1. //yy and cr form a smooth solution of (2.1|-(2.2), then the quantities 
defined by (5.1) satisfy the bounds 

Ait) < E2{t), (5.5) 

C{t) < E2{tfl^E:i{tYl\ (5.6) 

Similarly, suppose {r]i{t) , . . . , r]n{t)) and (cri(t), . . . , (7„(i)) form a solution of 
( |2.14 ) and (2.15) with the odd extensions (2.12). Suppose also that we have Ui = 
{V+77i, V+77i+i) ^ for \ ^ i ^ n — 1. Then a and c given by (5.3 1 satisfy the 
bounds 



a{t) < e2(t) 

c{t) < e,{tr/'e,{ty/\ 



(5.7) 
(5.8) 



Proof. We will just prove the discrete bounds for a and c in detail; the bounds for A 
and C can be proved using the same techniques, or we can view them as a limiting 
case of the bounds for a and c. 

The estimate (5.7 1 for a is easy: by Proposition 3.2 we have 



, n 



By Proposition 3.3 we have |(V_^iG)fej| ^ 1, and thus 



1 " 

(V_a),^-2|VH 



I'7ll5,l ^ 62 



by (4.22). We then have (5.7). The inequality (5.5) is proved identically. 

The bound (5.8) for c is more complicated. The first step is to differentiate the 
equation (2.15) in time to get 

-<V+ry,V_V+(aV+?7)> = <V+?7,V_V+((7V+r7)> + <V+?7,V_V+(f7V+r))> 

+ 2<V+77,V+fj> 
= 3<V+??, V_V+(aV+77)> + <V+??, V_V+(aV+77)>. 

Thus (T satisfies the same kind of equation as a with the same endpoint conditions, 
so we can use Proposition |3.2| to write 



1 ^ 

iV-&)k = - 2(V-aG)fe,(3<V+%,V_V+(c7V+?7),> + <V+r/,,V_V+((7V+?7),>). 
As above, the fact that |(V_^iG)fej| ^ 1 implies that c ^ |A|, where 

^=n^ (<V+77,,V-V+(aV+r,),> + <V+ry„V_V+(aV+77),>). (5.9) 



Applying the summation by parts formula (2.17) to this, using the endpoint 
conditions ?7„+i = r]n+i = fg, and performing some manipulations with the formula 
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<V+77, V+77> = 0, we get 

1 

^=n^ [v+'^,(3<V+%,V2+r7,>-<v2+r;,,V+%+i>) -4aJ<v2+r7,,v2+%>]■ 
i=l 

Using the bounds |V+(j| < 62 and < 62 from above, we obtain 
c ^ |A| ^ e2||?7||o,ilhllo,2 + e2||?7||i,2||'7||i,2 ^ ei\[(i2\[ez-, 



using Lemma 4.8 which gives (5.8). The proof of (5.6 1 is almost identical 



□ 



Lemma 5.2. Suppose 77 and a solve (2.1)-(2.2). Then B defined by (5.1) satisfies 
the estimate 

B(^)<^(l + ||r,||?>M?.^ (5.10) 

where Uo = JqI??*!^ ds is a constant depending on the initial condition. 

Similarly, suppose (7yi(t), . . . , ?7„(t)) and (cri(t), . . . , cr„(t)) form a solution of 



(2.14) and (2.15) with the odd extensions (2.12). Suppose also that |||?7|||3/2.2 ^ 



Then the quantity b in (5.3) satisfies the estimate 

=2|||r,|||l/2,2 



1 



b(t) < —e- 

Uq 



(5.11) 



where Uq is defined by (4.21) (and is constant due to Lemma J^.l). 



Proof. The bound (5.10) comes directly from (3.12): we have by Proposition 3.1 
that 



The fact that = 



Jo 



l'7t|lo,o 



< 



sx 

2 



x\ritxit,x)\'^ dx ^ 



i + hlll2 



The proof of (5.11) is similar, using (3.11) and (3.3) for the discrete Green 
function. The bound uq ^ WvWi 1 similarly follows from (4.12) since rin+i = 0. □ 



rjtWi I follows from (4.12), since rit{t, 1) = 



Example 5.3. In terms of the weighted energy (4.3), Lemmas 5.1 and 5.2 give 



upper bounds for A and B if the energy E2 is finite, while we only get an upper 
bound for C if E3 is finite. Simple examples show that these conditions are neces- 
sary: we can have Ei bounded while A and B are unbounded, and we can have E2 
bounded while C is unbounded. 

To obtain the examples for A and B, we consider (at time t = 0) the whip 
position 



?7(0,s) 



3(1 — s cos ( I Ins)) 3s sin ( I Ins) 



(5.12) 



/13 

which satisfies 7y(0, 1) = 0, |?7s| = 1, and \r]ss\ = This corresponds to a whip 
where the free end sits at ( —1= , 0) despite making infinitely many rotations around 
it (as s ^ 0). See Figur e [3} 

The Green function (3.2) can be computed explicitly to obtain a from \rjst\- If 
77(0, s) satisfies (5.12) and |77st(0, s)| = s^'^/^, one computes that Ei is finite while 
E2 and A are both infinite. If on the other hand \r]st{0, s)| = 1, we easily see that 
El is still finite while E2 and B are both infinite. 

The example for C is a bit more involved. Suppose 77 is given by (2.3), where 
9s{0,s) = s^^/^ and 9t{0,s) = s^^^^. It is easy to verify that E2 is finite at 
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Figure 3. The curve defined by (5.12), for which the curvature 
approaches infinity at the free end. Although the curve has length 
one, its free end wraps around the limiting point infinitely 
many times. We have plotted this heuristically in the inset, al- 
though the actual curve wraps itself up too tightly for these loops 
to be visible. 



this instant, while E3 is infinite. We have cr(s) = s at this instant by ( |2.2| , so that 
differentiating (2.2) with respect to time and using (2.1 ) gives atss{s) — s^'^/'^at{s) = 
— 3/s with boundary conditions (Tt(0) = and (Tsi(l) = 0. In this case we can verify 
that C is infinite. 

The fact that we cannot bound C unless is bounded is one of the main reasons 
why the energy estimates only close up at E^. Taking a time derivative of (4.1) 
as in (4.2) gives a number of terms of the form at which can only be bounded in 
terms of C, and thus in terms of E^^. 

We are now ready for an a priori estimate on the tension cr. Although the norms 
of a and CTs are easier to measure using the supremum, it is convenient to use 
weighted Sobolev norms for the higher derivatives of a. Thus we define the squared 
norm for a whip: 

m— 1 



1+^/2,1+2 



for TO ^ 1. 



£=0 



The discrete version is defined by the same formula: 

m— 1 

dm = 2 lkll£+3/2,£+2 for TO ^ 1. 
£=0 



(5.13) 



(5.14) 



Lemma 5.4. If a is a smooth solution of (2.2), then the norms (5.13) can be 
bounded by the energy (4.3) via 



Di < El, 
while for m > 5 we have 



Do < E: 



31 



and Dt, < E- 



3' 



Dm ^ Pm{Em-l)E„ 

where Pm depends only on E^-i- 



(5.15) 
(5.16) 



Similarly if a satisfies (2.18) with the condition ctq = 0, then the norms (5.14) 



can be bounded by the energy (4.18) via 
di < e\ 



.4 
-3' 



d. < e 



3' 



and dj! < e 



3' 



(5.17) 
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while for m > 3 we have 

where Pm depends only on Cm-i- 



(5.18) 



Proof. We will just prove the discrete estimates (5.17H(5.18); the estiraates (5.151 



(5.161 are proved using the exact same technique. The full proof is in Appendix 



A. 3 the basic idea is just to take iterated differences of (2.181 and estimate using 
Corollary |4.5| and Lemma 4.8 □ 



6. The main energy estimate 



In order to construct the solution of the partial differential equations (l.ll- 
1.2|, we want to find bounds on all the discrete energies (4.18 1-(4.19 ) which are 



independent of the initial conditions and of the number n of links. Then in Section 
[7]we will find a subsequence that converges to a solution. As a consequence, we can 
show that the motion of a chain converges to the motion of a whip as n approaches 
infinity, in the sense that position, velocity, and acceleration all converge. 
We now wa.nt to estima-tc the time evolution of the energy G'm- 

Our strategy 

will be to bound de^^/dt in terms of the energies e^; we will then use the fact that 
e„i and e,„ are equivalent (since a/s is bounded above and below by Lemmas 5.1 
and 5.2) to get an inequality for dem/dt in terms of Cm- In proving it we will use 



Lemmas 5.1-5.2 and|5.4| in an essential way. 



Theorem 6.1. Let n e N, and suppose (?7i(i), . . . , 77„(t)) and ((Ti(t), . . . , (T„(i)) 
form a solution of (2.14) and (2.15) with (7Q{t) = and rin+i{t) = 0, and that rj 



and a extend to sequences satisfying the oddness condition (|2.12 ). 



Then the energies (4.18) and (4.19) satisfy the estimates 



d£3 
dt 



(6.1) 



for some M3 independent of the initial data and ofn. In addition the higher energies 
satisfy 



dCyyi 
Ht 



< 



A^m(em-i)e^ 



(6.2) 



for every m > 3, where Mm depends only on Cr, 



Analogously, if rj and a form a smooth solution of (2.1) and (2.2), then the 



energies (4.1) and (4.3) satisfy the estimates 



dE:, 



(6.3) 



dt --'"^^ 

for some M3 independent of the initial data. In addition the higher energies satisfy 

dE„ 



dt 



< 



Mm{Em-l)En 



for every m > 3, where Mm depends only on E„i-i- 
Proof. The proof is in Appendix |A.4| 



(6.4) 



□ 



The fact that the energy estimates only close up at m = 3 is perhaps explained 
by the following observation, which is easier to understand in terms of the spherical 



representation (2.3) 
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Proposition 6.2. Let denote the unit hall in M"*. The pair {0,a) is a smooth 
solution of (2.4) -(2.5 I if and only if the functions ip: 13^ and a: D'^ — >■ M"'" 



2^ „..^_'1N!) 

4|x|2 



defined by 

(p{x) = d{\x\'') and a{x) 

are spherically symmetric solutions of the equations 

(fitt = div(agrad<^) 

(6.5) 

Aa- |grad(/9pa = -\(ptr. 

with Neumann boundary condition dyip = for if and Robin boundary condition 
dvOL + 2q! = for a on dD^ = S'^ . Furthermore any smooth solution has a > 
everywhere, so that the hyperbolic equation for ip is nondegenerate. 

Proof. Setting cr(s) = Asa{s) and changing variables by s = r^, we easily see that 



(12^-(12^ become 

Now the operator d'f + is familiar as the spherically symmetric Laplacian on 
W^, and hence we recognize both terms above as coming from the Laplacian on 
under the assumption that a and 6 are both spherically symmetric. The boundary 
conditions are easy to check. □ 

The fact that the degeneracy can be removed if we work in a higher-dimensional 
space, and thus in some sense the equations naturally "live" there, is essentially the 
reason why we need higher than usual Sobolev order for the estimates to close. 

7. Local existence and uniqueness of the solution 



Now we can finally prove the local existence theorem for the system (2.1 )-(2.2) 
of partial differential equations. The fact that Theorem |6.1| gives us estimates for 
Cm in terms of Cm that are independent of n allows us to construct the solution as 
a limit of a subsequence of discrete solutions as n — > oo , following the technique of 
Ladyzhenskaya [D and references therein. 

7.1. The discrete interpolation. We first need to establish the interior approx- 
imation of the space of whips by the space of chains, which allows us to go from 
estimates on e™ given by ( |4.18[ ) to estimates on Em given by (4.3) and back. 



Consider any function rj: [0,1] such that \ris\ = 1, with rj extending to 

an odd function through s = 1, such that the seminorms for 2 ^ £ ^ m are 

all finite. For each n e N we want to approximate 77 by a sequence rjk G K'' for 
1 ^ fc ^ n, extend it for fc > n by rjk = — ?72n+2-fe, have it satisfy |V+?7fc| = 1, 
and have uniform bounds on the discrete Sobolev seminorms Hf^H^.^ in terms of the 
smooth seminorms that are independent of n. 

The complication arises from handling the constraint \r]s\ = 1. Although it is 
relatively easy to approximate functions by sequences in the norms we need, the 
typical discrete approximation will not satisfy the condition jV+ryj = 1, which 
means it does not actually represent a chain. We deal with this by using the 
spherical representation r]s{s) = (cos9{s),sm9{s)) as in (2.3|. (Although this 



formula works only when d = 2, we can use a similar procedure in higher dimensions. 
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using generalized spherical coordinates.) Using 77(1) = 0, we can easily reconstruct 
?7 if is known. We can then approximate the function by a sequence 9k and 
rebuild rjk using the formula ?/„+i = and V+?7fc = (cos 0^, sin 0^) when d = 2, 
with a similar formula in higher dimensions. 

Fortunately, the Sobolev norms of r/ and 6 are closely related. 



Proposition 7.1. Ifrj: [0,1] 



elated to : [0, 1] 



by the formula (2.3 I, 



with 77(1) = ri"{l) and 0'{l) = 0, then houndedness of the squared norm 

^ = J' (s2|0'(5)|2 + s3|0"(5)|2 + s4|r(s)|2) ds 

is equivalent to houndedness of the squared norm 

B = J' (^s^\r,"{s)\^ + s^\v"'{x)\^ + 1 77(4) (x) 1 2) dx. 
Proof. We easily compute that 

h"(s)p = 0'(s)^ w"{s)\' = 9"{sr+0'{sr, 

and |?7(^)(s)|2 = {9"'{s) - 9'{sff + W'(sf9"{sf. 



(7.1) 



(7.2) 



(7.3) 



An integration by parts using 9'{1) = shows that (7.1) and \I.2\ are related by 
B = A + J' [s^9'{sf - 8s^9'{s)^ + 15s^9'{sf9"{sf + s^9'{sf^ ds. (7.4) 



Repeated use of the basic weighted Sobolev inequalities of Theorem |4.2| allows us 
to express every term on the right side in terms of A, so we get an inequality of the 
form B < A + + A^. 

In the other direction, (7.4 1 gives 

As;B + 8 f s^\?]"{s)\Us, 
Jo 

and again using Theorem 4.2 gives A < B + B^. 



□ 



We can derive the same sort of result for difference quotients. If we write V+Ty^ = 
(cos9k,sm9k) for 1 ^ A: ^ n, then the analogues of (7.3) are as follows: 

|V^,,P = 4n^sin^ (^), 

I V^f ^ 16n^ sin^ + ) sin^ (^) + sin^ (§) 

sin ( ^^^""ir-^^-^ ) - 3 sin (^i^) f 
+ 48 sin^ sin {^^^) sin (^i^±^^l?±^) . 



in" 



For sufficiently large n, we can proceed as in Proposition |7.1| to show that the 
discrete squared norms 

n— 1 n— 2 n— 3 



and 



n— 3 



-i ft — J. -J It — -J (t — 



k=l 



THE MOTION OF WHIPS AND CHAINS 



25 



can each be bounded m terms of the other. 

Thus for either whips or chains in two dimensions, we can work directly in terms 
of Sobolev norms of 9. The most convenient way to map from Sobolev spaces of 
continuous maps to Sobolev spaces of discrete sequences is to use orthogonal poly- 
nomials. (A direct approach, using values of the function on a discrete grid, does 
not work for our purposes since bounds on the differences require more smoothness 
of the function than we have.) 

The only complication is the oddness requirement on rj (and the discrete oddness 
criterion (2.12)). In terms of the spherical variable 6, oddness of 77 through s = 1 
translates into evenness of 9, i.e., there is an extension of 9 to [0,2] such that 
9{2 — s) = 9{s). Similarly the discrete oddness condition (2.12) translates into the 
discrete evenness condition 02n+i-fe = ^/c- These conditions are easy to handle if 
we extend the interval to [0, 2] (or extend the sequence to {1,2,..., 2n}) and use 
Sobolev seminorms with symmetric weights 

((9,9))pj = \ p{sy+^\9^^\s)\^ds where p{s) = s{2 - s) (7.5) 
Jo 

and 

«O,0\, = - 2 P^^^'^\KO,\' where = M^^ill^. (7.6) 

k=i 

These norms are clearly topologically equivalent to the weighted norms we have 
been using, where the weights are s and ^. 

Theorem 7.2. There are polynomials Qm{s) on [0,2] satisfying Qm{'^ ~ s) = 
Qm{s), and such that 

«Q.,Q™»P.=W™., ^here = _ ^^^^^^^^-^^^^ _ (7-7) 



in the weighted Sobolev seminorms (7.51 for all j ^ 0. Thus we can expand 9{s) = 
Yim^i^mQmis) and obtain 

m— 1 

for all j ^ 0. 

There are also, for each n e N, discrete polynomials qm{^) defined for 1 ^ fc ^ 2n 
and 1 ^ m ^ n, satisfying qm{ '^"\l^'' ) = 9m (^) md such that 

«'Zi?,'7m»pj = Simrmj, (7.9) 

where r,nj is as in (7.7 1. Hence if 9k = Xim=i '^mlmi^) for 1 ^ k n, then 

n 

«0,0X,, = J] r^.al. (7.10) 

rn—l 

Proof. The desired polynomials come from a slight variation on the classical Le- 
gendre polynomials and the Chebyshev polynomials of a discrete variable. (These 
are special cases of the Jacobi and Hahn polynomials respectively, with parameters 
a = /? = 0.) The desired formulas follow from general properties of continuous and 
discrete orthogonal polynomials; see Nikiforov et al. |NSU| for a good reference. 
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To obtain ([73, we set Q„,{s) = A„ P^,„_i(l-s), where P,(a:;) = 2^sF(a^'-l)" 
is the usual Legendre polynomial given by the Rodrigues formula and K^^ is a 
constant chosen to make Qm orthonormal when j = 0. To obtain (|7.9[), we set 



= ^[^^™-^i(^'2n+ 1) - ht"\{k-l,2n+l)], 

where 

h^^'°\x, N) = ^^(1 - E-^Y[{x + 1) • • • (a; + r){N - 1 - x) ■ ■ ■ {N - r - x)] 
r] 

is the Hahn polynomial given in terms of a discrete Rodrigues formula, with E 
denoting the integer shift operator, and again kmn is a constant chosen to give 
orthonormality when j = 0. 

Checking all the conditions is routine using the formulas in jNSUj. □ 

For each n 6 N, we can define the map Fn which takes a continuous angular 
function 9{s) to a discrete approximation 9k, and the map G„ which takes a discrete 
angular sequence 9k to a continuous angular function 9{s), by the formulas 



m— 1 m— 1 

n n 

9k = Y "-rnqmi^) '-^ ^ arnQmis) = 9{s), 
m — 1 rn — 1 

where the coefficients are obtained using orthonormality by 

Am= 9{s)Qmis)ds and a™ = - V 6'fcgm(f ) 



(7.11) 



By the formulas (7.8) and (7.10), we can bound the continuous and discrete Sobolev 
norms of any order j in terms of each other using this map. Furthermore G„ is an 
isometry, Fn o G„ is the identity, and G„ o Fn converges strongly to the identity as 
n ^ 00 in any weighted (p, j)-norm. 

Thus given an initial condition 77(0,5) = 7(5), we can write the discrete initial 
condition 7„ as 

(7«)fc = 2(V+7„)j(i) = {cos9k,sin9k) 

j=k j=k 



where 9k is the discretization obtained from (7.11). And conversely, if we solve the 
discrete chain equations to obtain rik{t), we can construct an approximate whip 
solution by finding, for each t, the angles 9k{t) and using ( |7.11 1 to obtain the 
function 9{t,s), then reconstructing rj(t,s) = ( cos ^?(i, x), sin 0(t, a;)) dx. 

We clearly have a similar construction for the velocity r]t{t,s) in terms of the 
angular velocity 9t{t, s), which works based on the formulas 

ilst{t,s) = {- sm9{t, s),cos9{t, s))9t{t, s) 

V+f]k{t) = {-siii9k{t),cos9k{t))9k{t). 

These constructions ensure that we can go back and forth between whips and 
chains while preserving the Sobolev norms as well as the constraint equation. 
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7.2. Uniform energy bounds. Now suppose that the mitial whip conditions 
77(0, s) = 7(5) and rjt{0,s) = w{s) have bounded energy £"3(0) given by (4.3), as 
well as satisfying the constraints |7'(s)p = 1 and (^j' (s) , w' (s)} = 0, and have odd 
extensions through s = 1. Using the procedure of the preceding section, we know 
that for each n e N there are discrete initial conditions 7„ and Wn such that the 
discrete energy 63 given by (4.18) is bounded uniformly, independently of n. These 
approximate conditions converge strongly in iV4[0, 1] and Af3[0,l] respectively to 
the actual initial conditions. 

Lemma 7.3. Suppose 7 and w are initial conditions as in Theorem \l.l\ and sup- 
pose discretizations 7„ and Wn are defined as in Section 7.1 Let (?7n)fc(i) and 
{o'n)k{t), for 1 ^ k ^ n, be the solution of equations (2.14) and ( |2.15 ) with 
VniO) = 7n and 77„(0) = w„. 

Then there is a T > such that the discrete energy e^lt) defined by (4.18) is 
bounded uniformly on [0,T] and uniformly in n. 

Proof. Since the discrete energy 63(0) is bounded uniformly for all n, we conclude 
by Lemma 4.8 that |||?7n(0)||| 3/2,2 is uniformly bounded for all n. In particular we 
know the hypotheses of Lemma |5.2| are satisfied for n sufficiently large. Fix such 
an n. 



By Lemma 5.2 we have 



Sk 



^b{t) 



< 



o2|ll'?(t)lll3/2 



Mo 



We want an estimate for the evolution of |||?7|||3/2 2- ^"^^ ^ ^ 
have by Proposition 4.4 and Lemma [4. 8| that 

fX''^"^"^\nVk{t)\') = 24^/^)<V^%(t), V^^%W> 

< llhWllli,2|||^( 

< 63 (i). 

Since this is true for any fc, we conclude 

L 



(7.12) 
1}, we 



2,2 



3/2,2 



3/2,2 + 4 



J 

Jo 



63 (r) dr 



(7.13) 



for some constant L, independent of t and n. This bound also ensures that the 
hypotheses of Lemma 5.2 are satisfied for sufficiently large n as long as e3(i) is 
bounded. 

By the definitions (4.18) and (4.19), we clearly have 



[l,b{t)}%{t), 



and we conclude by combining (7.12) and (7.13) that 



63 (f) ^ K exp 



63 (r) dr e3(t), 



(7.14) 



for some constant K which is also independent of t and n. 



Let y{t) = e3(T) dr and let z{t) = ^^{t) dr. Then (7.14) can be written as 

-^,„lt^dy ^ ,^dz 
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and integrating both sides yields 



-Lyit) 



^ 1 - KLz{t). 



(7.15) 



Now we use Theorem 



6.1 



to get ^ ^ M3e3(t)'' for some M3 independent of t 
and n. Using ( |7.14[ ) and ( |7.15 ), we have 

'dzV 
... dt) 

so we obtain 



^3^^^ ^^"^ {l-KLz{t)y 



dt^ ^ {l-KLz{t)y 
Dividing by (^)^ and integrating, we obtain 



z'{t) ^ z'(0) 



1 + J 1 



{l-KLz{t)Y 



-1/5 



(7.16) 



(7.17) 



where J 



6L 



Another integration gives a bound for z{t) on some time 



interval [0,T], which depends only on 63(0), and (7.17) gives a uniform bound on 
63 (i). Combining this with (7.14) and (7.15), we get a uniform bound on 63 (i) as 



well on the same time interval. 



□ 



Now having obtained a sequence of chain solutions r?„(t), bounded uniformly in 
the discrete weighted Sobolev norms uniformly on an interval [0,T], we use the 
technique of Section [TH] to interpolate. For each n we obtain an approximate whip 
solution ??„: [0,T] x [0, 1] M'* for which the energy E^(t) is bounded on [0,T] 
independently of n. We can then extract a subsequence which converges in the 
weak-* topology on L*([0, T], iV4[0, 1]). 

Before doing this, we prove one final lemma, a compactness result analogous to 
the usual Rcllich theorem. 



Lemma 7.4. Let A^m[0, 2] denote the space of functions 77; [0,2] 
the norm 



M 



JO 



ds 



such that 



(7.18) 



is finite 
Then N 



m + l 



[0, 2] is compact in A^„j[0, 2] for each m ^ 0. 



Proof. Expand 77(5) = ^j'^QWjPj{l — s), where Pj are the standard Legendre 
polynomials. Then as discussed in Section[7.1[ we have 



l^ll 



m r/j 

V 



I] 



Hence the embedding t: Nm+i — >■ N^, 
dimensional range, so it is compact. 



is a norm limit of operators with finite- 

□ 



As noted in Section [O] for functions on [0, 1] that are restrictions of odd func- 
tions on [0,2], the norm on iVm[0, 2] given by (7.18) is equivalent to the norm on 
Nm[0, 1] given by ( |1.3[ ), and thus we get compactness of Nm+i[0, 1] in A'^,„[0, 1] for 
functions with an odd extension through s = 1. 

We now establish the existence part of Theorem |1.1| 
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Theorem 7.5. Given initial conditions 7 and w as in Theorem there is a 
T > such that there is a solution rj of the system (7.19) in L'^^(\Q,T\, N^\Q,1\) n 

w^i'*([o,r],7V3[o,i]). 

Proof. For each fixed t and each n e N, construct a continuous approximation of 
the chain 77,1 (i) as in Section |7.1[ and call it Then by Lemma 7.3 we get 

a uniform bound on on E^{t) in some short time interval [0,r]; in other words, 
the family ^„ is bounded in L*([0, T], 7V4[0, 1]) n T^i'*([0, T], A^3[0, 1]). By the 
Alaoglu theorem, there is a subsequence that converges in the weak-* topology 
to 7? 6 L*([0,T],7V4[0, 1]) n VFi'«([0,T],iV3[0, 1]). 

By the compactness Lemma 7.4 there is a sub-subsequence 77^^. . which converges 

strongly to 77 in L*([0, T], iValO;!]) W"i'*([0, T], TVsLO, 1]). For any e > the 
convergence is strong in H^\e, 1], and thus by the usual Sobolev embedding theorem 
also in C^[e, 1]. So we can take the limit of the system (2.14) and (2.15) pointwise 
to see that we have a solution of ( |7.19 1. □ 

The fact that all the estimates close up at the level of i?3, with all other energies 
satisfying linear differential inequalities, implies that the only way a solution which 
is initially can fail to be G^- for all time is if S3 becomes infinite in finite time. 
This gives a crude blowup criterion. 

Corollary 7.6. Suppose rj, a is a solution of the system 

Vtt = ds{aT]s), - Ivssf = -\Vstf , IVsl'^ = ^, 

7;(i, 1) = 0, (7,(1) = 0, c7(0) = 0, ?7(0,s)=7(s), 77t(0,s) = w;(s), (7.19) 

where we assume that 7 and w satisfy the conditions of Theorem \l.l\ 

Assume that in some time interval [0, T], the energy E^^it) is hounded uniformly. 
Assume further that Em{0) is bounded for all m > 3. Then Emit) is also bounded 
in [0, T] for all m > 3. 



Proof. By equation (6.4), we have for k ^ i that 

dEk 



dt 



< 



Mk{Ek-i)Ek. 



Furthermore since E2{t) is bounded, so is B{t) 
thus 

dEk 



sup 



s u{t,s) 



by Lemma 



5.2 



and 



dt 



< Mk{Ek-i)Ek 



for some function M^. So by Gronwall's inequality, Ek{t) is bounded on [0,T] in 
terms of -Efc(O). Thus finally Ek{t) is also bounded on [0, T]. □ 

We now complete the proof of Theorem |1.1| by proving uniqueness. 

Theorem 7.7. Suppose 7 and w are functions on [0, 1] as in Theorem If 
(771, CTi) and (772,(72) are two solutions of (7.19), both in 

L*([0,r],iV4[0,l]) n W^i-«([0,r],7V3[0,l]), 

with the same initial conditions 

771(0, s) = 7/2(0, s) = 7(s) and dtrii{0, s) = dtri2{0, s) = w{s), 

then 771 (t, s) = 772 (t, s) and (7i(t, s) = (72 (i, s) for all t 6 [0, T] and all s 6 [0, 1] . 
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Proof. The proof relies on an energy estimate for the differences at the level of the 

first energy, 

Eibli - = J (\dtm - 8tm? + s\d,'qx - 5,772p 

+ s\dtdsm - dtdsm? + s^\8lm - 5^??2|^) ds. 



We estimate this energy using a Gronwall inequality, as in Theorem 6.1 and Lemma 



7.3 The reason this works is that since r]i — r]2 satisfies a linear PDE whose 
coefficients involve the known quantities 771, 772, tri, and <72, we can use Corollary 
4.5| in a more effective way to put all the weights on the known terms. The full 



proof appears in Appendix |A.5[ □ 
Finally we discuss some refinements of these results. First, given a solution 77 



of (1.1 1 with finite, we can check using the differential equation that |j5(?7||2,2j 
ll^^jylli 1, and ||5j^77||o.o can all be bounded in terms of E^. Hence the solution is 
also in V2'«([0, T], iV2[0, 1]) n M^3'*([0, T], iVi[0, 1]) n VF4,oo([o, T], iVo[0, 1]). Now 
a well-known general technique (see e.g.. Lemma 11.9 of jRRp shows that rj is 
continuous as a curve in A^4[0, 1], 77 is as a curve in A*'3[0, 1], etc. 

8. General remarks and future research 

In this paper we considered the whip with one fixed and one free end as boundary 
conditions. The other possibilities are to have two free ends, to have two fixed 
ends, and to have periodicity. All of the estimates in this paper have analogues 
in those cases. When there are two free ends, the tension must satisfy cr(0) = 
and cr(2) = 0, so the appropriate weighted norms look like the square root of 

5*^(2— s)*^ l/^'^) (s) P ds. Since we have essentially solved the problem with one fixed 
end by constructing an odd extension in order to turn the problem into a string 
with two free ends on [0,2], we expect that the same estimates prove existence 
for an inextensible string with two free ends. When there are two fixed ends, or 
when the whip is periodic, the problem becomes simpler since we can use ordinary 
Sobolev spaces for the estimates. In this case we expect the energy estimates to 
close up at the level of 62 rather than 63. 

The addition of gravity brings some complications. One is that the boundary 
conditions change, and oddness through the fixed point is no longer enough to 
satisfy the conditions automatically. (This is already an issue even for the wave 
equation with constant coefficients, if an external force is imposed which does not 
respect the boundary conditions.) The other complication is that if the whip is 
above the fixed point, the tension may become negative: the effect of gravity is 



to change the boundary condition in (2.2 1 to as{t,l) = (g,r]s{t,l)}, where g is 



the gravitational acceleration vector, and if as{t, 1) < then it is possible to have 
a{t, s) < for some t and s. In that case the evolution equation becomes elliptic, 
so the discussion becomes much more complicated. 

The blowup criterion Corollary |7.6[ that a smooth solution remains smooth up 
to time T iff supp^^^y i?3(t) < 00, can certainly be improved. Once we know a 
solution exists, we can use alternative methods to get better a priori bounds on it. 
Thess et al. have speculated that blowup for the periodic loop might be controlled 
by the U^' norms of \r]ss\ and |77st|, analogous to the way blowup for the ideal Euler 
equations is controlled by the L'^~' norm of vorticity. This is an interesting problem 
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to study, since we have a much greater handle on all aspects of this one-dimensional 
problem. We will explore this in a future paper. 

In addition, the geometry of the space of inextensible curves is interesting in its 
own right. Although the geometric objects are not smooth in the Sobolev topology 
(unlike on the group of volumorphisms) , the curvature formulas still make sense, 
and one can compute formally that all sectional curvatures are nonnegative. We 
can thus try to study stability of the motion from the geometric point of view (as 
in |AKj ). as well as the geometry of blowup. See [P] for details on this. 

A similar problem in higher dimensions is given by the motion of a flag attached 
to a pole in 3-space. Here our configuration space would be the space of maps of 
a rectangle into M.^ which are isometric immersions with one side of the rectangle 
held fixed. We expect to see a similar nonlocal coupled degenerate system, the only 



obvious difference being that the ordinary differential equation (1.2) becomes an 
elliptic equation in the spatial variables. 

The whip-chain equations are interesting partly in and of themselves, but espe- 
cially as a "toy model" of inviscid, incompressible fluids. There are some structural 



similarities between the equations (l.ll and (1.2 1 and the Euler equation for an 
ideal fluid, given in Lagrangian form by 

7jtt{t, x) = - gradp(i, T]{t, x)) and Ap = -TT{[Df]t{t, x) o r]^^{t, x)f) , 

with some boundary condition to determine gradp uniquely. Both systems involve 
a hyperbolic evolution equation for a constrained function, where the right side is 
given in terms of a function determined by a purely spatial differential equation. 
The technique of approximating a continuous system with a discrete system pre- 
serving the geometry may be interesting to apply to fluids directly. For example, 
in two dimensions we could consider a rectangular grid on a torus, the vertices of 
which are free to move as long as all quadrilateral areas are preserved. Although 
such a model may not have global existence (as edges of a quadrilateral may col- 
lapse to give a triangle without changing the area), we might still get some useful 
insight out of it. 

Appendix A. Longer proofs 
A.l. Proof of Proposition [375| 

Proposition. Suppose Gkj, r/k, otk, and j3k are defined as in Proposition 3.2 

- frt r cn irrt i ^ C ( r\ - 

^|V^77feP^u /or lsCfcsCn-1. (A.l) 



Assume the rjk are such that, for some v 6 (0, ^^], we have 



Then for every 1 ^ j,k ^ n, we have 



jk 



^ e"^^. (A.2) 



If G solves (3.2) and rjss is a smooth function, then we have 



r 

Jo 



inf ^ JLL ^here g = | s|?7s.f ds = MU,^. (A.3) 

o^s,x^i sx 1 + g 



Proof. Our strategy for proving (A.3) will be to first show that the minimum 
^ and ^ 

kj sx 



of the ratios ^ F-^^ and ^^^'^^ is attained at the ofF-diagonal corners; that is, 
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minis;,- /js;. 



nGin and info<s 



G{s,x) 



lim 



s->0 



G,(0,1). The 



proofs are nearly identical in both cases, so we will just give the discrete proof. 
Then we estimate the size of this value; here the proofs are different, and we can 
get a sharper estimate for the continuous case. 

2 

We first define a matrix F by F^j = j^Gkj- Clearly F is symmetric since G is. 
We want to prove that Fkj ^ Fi„. Note that for 1 < fc ^ n we have 



Fkj — Fk-i. 



[{k-l){Gkj-Gk-i,j)-Gu-i.j'\. 



(A.4) 



jk{k-\) 

First we show that we can decrease F^j by increasing the larger index, li k > j 
we know from Proposition 3.3 that (V_.iG')fcj ^ 0, so that Gkj ^ Gk-ij, and thus 
Fkj - Fk-i,j s; 0. Thus 

Fkj ^ F„j if fc ^ j. (A.5) 
Next we show that we can decrease Fkj by decreasing the smaller index, which 



is a bit more involved. Inspired by ( |A.4[ ), we define for 1 ^ k ^ n the auxiliary 
quantity Hkj = (fc — l){Gkj — Gk-i.j) — Gk-ij', then it is easy to compute that 
{V+^iH)kj = fc(V_4V+,iG)fcj, and we conclude using (3.8 1 that if fc < j then 
Hk+i.j — Hkj ^ 0. Since Hij = 0, this shows that Hkj ^ as long as fc ^ j. 
Then since Fkj - Fk-i,j = jk(k-i) Hkj for fc > 1, we have Fkj ^ Fk-i,j as long as 
1 < fc ^ j, and hence 

Fkj ^ Fij for 1 sc fc ^ j. (A.6) 



Combining (A.5) and (A.6), and using the fact that Fkj = Fjk, we obtain 

Fm. (A.7) 



min Fkj 

l^j.k^n 



We finally want to bound Fin from below. Using the formula (3.4 1 we have that 

ji-i 

r)i 1 r)i I 

Fin 



PllPln 



^ 11 — 1 



/3l Pi (3m+ 



(A., 



as 



(A.9) 



It is easier to estimate sums than products, so we rewrite (A. 

n—l n—1 

In Fin = 2 Inoffc - 2 ln/3fe, 
fe=i fe=i 

recalling that /3„ = 1. 

First we get an upper estimate for Xife^il^/^fc- Rearranging (3.5) and using 
ak = l- 2^|V+?7fep, we have 

\^lvk\' mvk\' 



Pk — Pk + l 



k+i 



1) + 



4n4 



Recalling that 1 ^ (3k ^ 2 for each fc, we conclude 



1 



pk-Pk+i ^ :%|v2+77fcp, 



and since = I7 we find 



n—l n—l 

Pn+J] iPk - Pk+i) 1 + ^ 1] iV'+r/fcl 



(A.IO) 



k^j 



k^j 
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Now since /3j ^ 1, we have ln/3j ^ /3j — 1, so that (incorporating the assumption 



2 ln/3fc 



n-1 



n-1 



2 |2 



2 |2 



k = l 



s; u 2 + 



C(l/2) \ 
\/n J 



(A.ll) 



where C(l/2) * —1.46 is a value of the Riemann zeta functionj^ 
Next we get a lower estimate for Xil-^il'^Q^fe- Since = 1 



2„2 , i:ne 
2ni/2fc3/2 ■ Since we have assumed u ^ 



assumption (A.l) yields ^ 1 
have afc ^ 0.8 > for all 1 ^ fc ^ n — 1. Now we want to get a lower bound for 
In ak ; this is a bit more delicate than an upper bound for the logarithm of In (3^ ■ 
Define c = -2C(5)/C(i) « 1.118 in ter ms of the Riemann zeta function. It is not 
difficult to verify that when 0.8 ^ ^ 1, then Ina^ ^ — c(l — Qffe). Thus we have 
Inafc ^ -(cu)/(2nV2fc3/2) for every fc, from which we conclude 



n-1 



n—1 1 



2nV2 ^ fc3/2 

k—l 



2nV2 ^fn ■ 



(A.12) 



Combining ( |A.11[ ) with (A.12) and plugging into (A.9|, we obtain lnFi„ ^ —2v. 
Using (A. 7), we obtain (A. 2) as desired. 



Now we will just sketch the proof of (|A.3|). We similarly establish that the 

;h woi 



G(s x) ' ' 

infimum of ^ ' ' is attained when s = and x = 1, which works the same way as 



Gs(0, 1). Letting 
(A.13) 



1. 



in the discrete case. So we just need to estimate lims^o 
J(s) = G(s, 1), we see that J satisfies 

J"{s)-W'{s)\'j{s)=0, J(0) = 0, J'(l) 

The minimum is then Gs{0, 1) = J'(0). 

Set A(s) = ln[J(s)/s]; then we can verify by explicit computation that (A.13) 
can be rewritten in two ways: 

d 



ds 



[s^X'(s)] = -s'^X'isf + s'|?7"(s)P, and 



ds 



[s{l - s)A'(s)] + A'(s) = -s(l - s)A'(s)2 + s{l - s)\r]"{s)\' 



Integrating the first equation from s = to s = 1 gives 

X'{1)^ C s'W'{s)\^ds^g, 
Jo 

and integrating the second from s = to s = 1 gives 



A(l)-A(0)iC r s(l-s)|77"(s)pds^ 
Jo 



(A.14) 



(A.15) 



Since A'(l) 

combine to give (|A.3[). 



J(iV 



■1 and A(1)-A(0) = In [J(1)/J(0)], estimates ( |A.14| and ( |A.15[ ) 



'^Having a precise estimate of this remainder is useful to make part of ( |A.11| cancel out | |A.12| , 
in order to make the estimate ||A.2| independent of n and thus a bit more elegant. 
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A.2. Proof of Theorem [13 



Theorem. Let f: [0, 1] R'' be . Then for any r > the norms (|4.4|) and 
(4.5) satisfy the weighted inequalities 

ll/llr-l.m ^ ll/llr,m + ||/|lr+l,m+l and 



< 



11/11™ + 11/11 



r+l,m+l • 



(A.16) 
(A.17) 



If in addition we have = 0, then these inequalities can be simplified to 



ll/ll'-i,™^ I 



/llr+l,m+ 



i and 



(A.18) 

nr,m ~ ll/llr+l,m + l- (A. 19) 

// / is instead a sequence {/i,--- ,/«} with values in M'^, then the inequalities 



(A.16)-( A.17) also hold if the norms are interpreted as (4.7) and (4.8), while the 

0. 



(m) 



inequalities (A.18)-( A.19) hold if f!^_ 

Proof. It is clearly sufficient to prove these inequalities when m = 0. To derive the 
discrete versions, we do the following. Let p be any real number. Then for any 
fc 6 {1, . . . , n — 1}, we have 



\{p + k)fk+i 
so that 



k{k+p) 



V+/fcP+pfc(|/fc+ip-|/fcP), 



/fc+il -l/fep). (A.20) 
Note that if we define /o in any way at all, the equation is still satisfied at fc = 0. 



Furthermore we have the easy-to-verify formulas ks 



(9-1) 



(9-1) 
k + l 



ns[^'^ and V+s^'^ 



which are valid for ^ 0. 



For any real q > 0, multiply (A.20) through by s^i^_^^^ and simplify to get 



Now notice that the last term simplifies to 

nps['\\f,^,\'-\f,\')=p[wM'^\f,r)+n{si'^ 



„(9) 



l)l/fe+lP] 



using V+s^''^ 



qs^^_^_P . Thus we have 



Now let i and j be any integers with ^ i < j ^ n. Summing all the terms from 
k = i to k = j — 1 and using the telescope formula ^ 21=1 V+fefc = bj — bi for any 
sequence {6^}, we obtain 



^ 



pip - q) 



49-l)|^^|.^lgfc±P.(9), 

A; — i 



+p{s^;^\h? 



l/d'), (A.21) 

after reindexing the first sum on the right side. This is the basic building block for 
all the other inequalities in this proof. Now we consider some special cases which 
will together prove (A.16)-(A.17). Take any r > 0. 
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For any integer i with 1 ^ i ^ n, if we set p = q = r and j = n in (A. 21 ), 
we get 

1 h _i_ 

k — i 

and since ^^s^J'^ = s^^'^^\ we have 



(A.22) 



We use this to obtain (A. 17); if /„ = we obtain (A.19I. 
Next, if we set p = |, g = r, i = 0, and j = n, then we get 



fc=2 fc = l 

Noting that fc + | < fc + r for any r > and any k, we obtain after solving 
for that 



<: 



r-1,0 ^ ^2 II-' llr+l,- ^ 



.l + '^'-Vnl 



(A.23) 



which is used to bound (A. 16 1. If /„ = we obtain (A. 18). 

Finally we get an upper bound for |/„|. Choose q = r +1 and p = — (r + 1) 

with z = and j = n. Then we have 



k — r — 1 



^^2(^ 1 

n f—', n n 



Xr+l) 



fe=l 



from which we conclude 

4'-+i)|/„P^2(r+l)||/||2„+-A_||/||2^^_^. 

Now we obviously have 

"+'"1^112 



(A.24) 



lr+2,1 



■l/ll 



r+l,l5 



and plugging into ( A.24[ ) gives 

+ 4r + 1 



11/11^+1^1 + 4(r + 1)11/11^,0, 



(A.25) 



r(r + 1) 

Combining ( |A.23[ ) with ( |A.25[ ), we g et ( |A.16p . Combining ( |A.22p with ( |A.25p 
and taking the maximum, we get (A.17). □ 

A. 3. Proof of Lemma 15.41 



Lemma. If a satisfies (2.181 with the condition ctq = 0, then the norms (5.14) can 
be bounded by the energy (4.18) via 

di < e|, d2 < e|, and d^ < e^, 

while for to > 3 we have 

where Pm depends only on Cm-i- 



(A.26) 
(A.27) 
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Proof. Applying the shift operator E to (2.18), we obtain 



Via 



;,|2 



Thus as a first step, we have 

Ik IU+3/2/+2 = ||V + Cr||£+3/2/ 

= \lE'a\EWlr^\' + la\Vlrf - {EW+^f 

< ikiv^7,ni,+3/2^, + iiiv+77ni,+3/2, 



^+3/2/ 



using the fact that is a bounded operator in any norm. (The technique we use 
will make it clear that the norm of E'^a\E'V'^r]\'^ is comparable to that of crlV^ryl, 
so there is no reason to study it separately.) 

To simplify notation a bit, let / = jV^Typ and g = |V+77p. Then the inequality 
above is 



l'''ll£+3/2,f+2 ~ 



k/ll 



^+3/2, 



+ \\9\ 



e+3/2,i- 



(A.28) 



Our first goal is to bound ||o-/| 



e+3/2,e 



in terms of the norms of a and the norms of 



/. Note that we have bounds on a and V+u in the maximum norm by Lemma 5.1 



while for higher differences of a the bounds are expressed in terms of Euclidean- type 
norms. So the complication comes from taking this into account. 
Using the general product formula (2.16) for differences, we have 



£+3/2,e 



1 " 

^ k = l 
1 " 



k = l 



e ^ n 



J=0 



Now the exceptional cases are when j = £orj = £— 1, because there we want to 
use the maximum norm on a directly. In the other cases, we still need to use the 
maximum norm, but we will give it some extra weighting so that we can use (4.11 ). 
So we have, using — ^ a and |V+it| ^ a, that 



k/ll 



< 



.,2 n 



,(f+3/2) 2 



k+el'^ + fkl + 



+ > ( max s 



(j + 2) 
k 



n (e+3/2) 



(1 " 



(i+2) 
fc 



Using the bounds s^^'^'^^^^ s1_^_g 



4^+^/^) and 



\j+2) ^ '^k+j from Proposition 



4.4 we obtain 



£+3/2/ 



^«^ii/ip 



i+7/2,e 



+ «^ll/ll^+3/2,£-l + Xi 
J=0 



2 I 

J+2jl 



-3-1/2,1- j- 



(A.29) 
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A straightforward apphcation of (2.16) and the basic estimates of Theorem 4.2 
proves the inequalities 



\l+7l2,i ^ eae^+i + for ^ ^ 0; 



l/ll£+3/2,£-i ^ eae^+i + ej for l^l; 

< eaCj+a + e|+2 for j ^ 0;. 



£+3/2,£- 



To conclude, we need to estimate the other term in (A.28). We will show 

|2 



(A.30) 
(A.31) 
(A.32) 

ll5ll'+3/2/ ^ eae.+ i + e,^ for i ^ 0. (A.33) 

We will see that the extra half-power in the weighting is only necessary for the 
norm of g; otherwise we could have worked with ||/||^_|_3 etc. instead. We have 

ll5ll,%3/2,. ^ 2 ■.„ll<Vi+^),i?^ Vt+l-^'^>||,2^3/,,o 

~ lll'?lll 1/2,1 + = l lll'?llli + lj + lll'?ll£-j,«+l-j 

< 636^+1 + Xijil ej+2et+2-] ^ 636^+1 + ef. 
Plugging ( |A.30P -( |A33| into ( |A.29[ ) and ( |A.28| , we obtain 

lkllf+3/2,^+2 ^ (l + a^)(e3e£+i + e^) +2^10(636^+3 + e2_^2)lkll 



-j-l/2,^-j- 



(A.34) 



Recall from Lemma |5.l| that a < 62 , so that 1 + < e| . So when ^ = or 



^ = 1 we get ||cr||3y2 2 ^ 6^6361 ^ e| and 110-115/2 3 ~ 636362 ^ e\. From this we 
have c?i < e\ and 1^2 ^ 6|. Now we can derive from (A.34) for m ^ 3 the recursive 
inequality 

dm ^ 62636m + e2e^j_i + 2^1^^(636^ + 3 + 6^+2)c'm-2-i- 

Plugging in m = 3 gives the base case ( |A.26 ), and induction on m gives (A. 27). □ 
A.4. Proof of Theorem ETTl 

Theorem. Let n e N, and suppose (771 (t), . . . ,rin{t)) and (cri(i), . . . ,CT„(t)) form a 
solution of (2.14) and (2.15) with (To(i) = and rin+i{t) = 0, along with the odd 
extensions (2.12). 



Then the energies (4.18) and (4.19) satisfy the estimates 



d£3 

dt 



(A.35) 



for some M3 independent of the initial data and ofn. In addition the higher energies 
satisfy 

^ < Mraiera-l)em (A.36) 

for every m > 3, where Mm depends only on em-i- 



Proof. As with the proof of Lemma |5.4[ the estimates for the whip and chain are 
proved in the same way, so we will just focus on the harder case of the chain 
(where nontrivial technical issues such as Lemma A.l arise). The essential step is 
the discrete analogue of the computation (4.2), together with the integration by 
parts employed to cancel out the highest-order term. Then we simply estimate the 
remainder terms using Corollary |4.5| and Lemma |4.8| 
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The first step is just to differentiate. We deal with, the terms in (4.191 one at a 
time. So fix an integer £ ^ 0. Then 



d 1 
dt n 



E (-f|Vt'7.|V.r^)|Vr^.|^) = ^ E (I + 2II) (A.37) 



where 



and 



dt 



dt 



(A.38) 



(A.39) 



k 

dt 



For ( |A.38[ ), if £ = 1 then we have ^ ^ csk, while if f ^ 2 the derivative "^"^ 
is a derivative of a product of i terms. Using ^ as*: and ^ cs^j, we clearly 
have 



dt'''^ 



^ £a'-'cs 



which obviously holds in addition for ^ = and £ = 1. Thus the sum of (A.38 1 is 
bounded by 



k=l 



^ a'-'cMli + a'cM^^,^,^, < e'+'e, for £ ^ (A.40) 



using Lemma [5T 



So our primary concern is (A.39). By analogy with the technique used to 
derive (4.2), we want to pull out the worst terms of (A.39) and collect them 
into a single exact difference which will sum to zero. It is easy to check that if 
i/j = cr(^+i)<£:V^?7, V++^r;>, then 

V_V = (ii;-V(^+i))<Vtr), V^++iV_ry> 

The middle term of V_V is precisely the second term of ( |A.39[ ), so we want to 
show that the difference of the remaining terms is relatively simple]^ We will do 
this computation in Lemma ] A. 1[ 

To check that this all works, we observe that the backward difference V-V' sums 
to zero: we have 

" k^l 



^The complication in this computation is the fact that rj — V-((TV+r;) — B"^ V+((TV+r)): 
we have nothing but V-i- in the rest of the formula, so the appearance of one V_ operator (or, 
equivalently, of one backward shift E~^) necessitates rederiving the formulas to get rid of it, rather 
than using a formula like ( |2.16[ l directly. The reason we don't want to see a mix of V_ and V-i- 
operators is because later in the proof we will need to use the fact that {V+t;, V+J?) = 0, while 
there is no simple formula for {V_?;, V-i-?)). 
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which can be checked to vanish due to ctq = and the oddness conditions (2.12) 



(Recall that this is precisely the reason that the summands in our energies (4.181 
(4.191 all terminate at n — [i/2\.) Hence we have 



.-[£/2j n-[i/2\ 

2 (n)fc= 2 (n-v-^)fe. 



Now by Lemma I A. 1[ we have 



(A.41) 



where the remainder terms are given by 



3=0 



Therefore if £ ^ 1 we have 
^-[e/2\ 



' E (11). 



k=l 



n-[e/2\ g 



where 



e+i- 



(A.42) 



(A.43) 



(A.44) 



for I ^ i ^ £. Here we use the fact that by Proposition |4.4[ the shift operators 
E^^^ that appear in (A.42 1 are bounded for any j; since having or not having the 



shift operators attached to a doesn't change the estimates in any way, we may as 
well ignore them. (Having the shift operators attached to one of the 77 terms would 



cause a problem, which is why we need Lemma A.l 
Obviously when 



1, we must have i 



case every summand in (A.44) involves (V+r/^, V+?7fe) = 0, so that 5*1 



1 also in the sum (|A.43|), and in this 

if £ 



1. 



Hence we will assume £ ^ 2 to estimate Si 
We will show that for any 1 ^ i ^ £, 




if ^ = 2 or ^ 
if ^ ^ 4, 



(A.45) 



for some function Li^g of e^-i. The bounds for 1 ^ z ^ £ — 1 are all basically the 
same, while the bound for the i = i term requires another trick. 



Using Corollary |4.5[ and Lemma |4.8[ it is straightforward to verify that S^^ 

for 2 ^ i ^ £ — 1. Then using Lemma 



< 



Vdse^ and that 5*4 < ^ee+2-i\/di+i 
we obtain (A.45 ) when 1 ^ i ^ £ — 1 



5.4 



The last case in (|A.45|) is when i 



£. (Recall that we are assuming £ ^ 2 since 
1, we have {V+r/, V+?)) = 0. Applying the 



5i = when £ = 1.) Since |V. 

difference operator V^^^ to both sides and using the product formula (2.16), we 
get 



p=l ^ 



p 



(A.46) 



40 



STEPHEN C. PRESTON 



Thus 



£-1 n-{t/2\ 



, n 
=1 



We obviously want to use the Cauchy-Schwarz inequahty on this, and to do so 
we have to assign the weight s). to the individual pieces. Recall that we have an 
estimate for ||cr||£+i/2,f+i ^ \A5 ^ ^J Piie-e-ij^e from Lemma 



5.4 



which means we 

must pull the power 3(^+^/2) out with this term, leaving us with s^^^^/^) for what 
remains. The worst term is the one with p = £ — 1, for then we have to estimate 
V+77I • I V+'^yl ||£_i/2 Q- want to pull out the supremum norm of V+?7, but we 



need some positive weight on it in order to be able to use ( 4.11[ )[^ 
Using Corollary |4.5| again, we compute 



5.<i];:;iii?-iv^/vii,+i/2,o 



ii?fvt-^77iivr\i 



£-1/2,0 



< Ik 



£+1/2, £+1 



E'-'y+v\\\i/2.4yiv\U-ifi 

+ 2^:1 ll|V^/S|||,,o||i?^vt''^||.-p- 1/2,0^ 

^ V3£(|||?7||| 1/2,1 hlU-l,£ + I]pIi lll'7lllp,p+l||'?ll^-p-l,f-p) 

We easily estimate the quantities here, using Lemma [4. 8[ to get 



p=l y 



Putting this together with the bound for from Lemma 5.4 we get (A. 45 1 for the 
cases i = I. 



Now plugging (A. 45 1 into (A.43), we get 



E (11). 

fc=i 



< 



< < 



'0 


e = i, 






63 


e = 3, 



(A.47) 



.i^£(e^-i)e£ ^^4, 



for some function Kg of ei- 



Using (A.40) and (A.47) in (A.37), we get 



d 1 

dt n 



k = l 



(^+1)1.7^+1^^12 



< 



4 ^ = 0, 

ei i=l, 

el £ = 2, 

el £ = 3, 

Me{ee-i)ee I ^ 4. 



Now summing from £ = to £ = m, we get (A. 35 1 and (A. 36 1 



□ 



To complete the proof, let us establish the formula (A. 41) 



^This is the only place in the paper where we actually need to split the weight into noninteger 
powers to make the estimates work. Without doing this, w e ca nnot close the estimates at the 
level of 63. This is why the Sobolev norms of u from Lemma 5.4 are defined the way they are. 
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Lemma A.l. We have the formula 



where Ri is given by (A.42). 



Proof. The product formula (2.16) yields 



Then using = 1 + ^ and some binomial expansions and identities, we obtain 



(A.49) 

We then compute the inner product of all terms with V^?). 
In this last sum of (A.49), notice that when i = we get 

using the obvious telescoping, which finally yields (A. 48 1. □ 
A.5. Proof of Theorem [73 

Theorem. Suppose 7 andw are functions on [0, 1] as in Lemma \7. 3[ //(771, cti) and 
{r]2,(J2) are two solutions of flA^ in L*([0, T], iV4[0, 1]) n VKi'*([0, T], iV3[0, 1]), 
with the same initial conditions 

771(0, s) = 772(0, s) = 7(5) and ^4771(0, s) = ^(772(0, s) = w(s), 

then 771 (t, s) = 772 (t, s) and (Ti{t, s) = (J2{t, s) for all t 6 [0, T] otic? all s 6 [0, 1] . 

Proof. Define the differences by e = ^(772 — 771) and 6 = 5(172 — fi), and the averages 
by = 5 ('71 + '72) and a = \{(Ji + (T2). It is easy to compute that these quantities 
satisfy the equations 

6ss = {\riss? + \£ss?)5 + 2<e,„^,,>?f- 2<e,t,^,t> (A.50) 
ett = ds (ctEs) + ds {6ri^) . (A.51) 

The quantities a and S have the same boundary conditions as those for cti and 172; 
similarly f] and e must both be odd through s = 1 since rji and 772 are. Furthermore, 
the fact that |5s77ip = 1 = |5s772 p implies that (es,r]g} = 0. 

We now estimate the norms of these quantities; the primary goal is to estimate 
the norm of e, but we will need the norms of the other terms to do this. For 
this purpose, we generalize the quantities A, B, C in|5.1[ the quantities from 



(5.14), the quantities Em from (4.3), and the quantities Em from (|4.l|): we will 

Em[e\ 



denote A[a\ = supo<;s<;i ^^v^, Em[e] = TiT=o 
time-dependent energy Em we use a for the weighting. 



£t\\h_ + ||e||^+i,^+i, etc. For the 
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It is easy to verify that 

A[a-\ ^ \A[a{\ + \A[g21 
C[(f] ^ \C[a{\ + \C[a2\, 
E^[7f\ s; i;,„ [771] ■ 



Bp] sc B[cfi\ + B[a2\ 
D^p] sc D,n{(T{\ + D 



We could proceed by imitating the proof of Theorem 6.1 to get a bomid for 
the energy £'2[£]; however it's simpler to use some alternative techniques to get a 
bound for i?i[e]. The reason this works is that we can separate all the estimates 
into low-derivative norms of e by compensating with high-derivative norms of rj. 

First we note that since 5{t^ 0) = 0, we have by the Cau chy- Schwarz inequality 
that (5(t, s)^ ^ s 5s{t, s)^ ds. In the notation of Definition 4.1 we conclude that 

(A.52) 



11-1,0 s; ll-^llo,!- 

Since <5(t,0) = and (5s(i, 1) = 0, we have using ( |A.50 ) and Corollary 4.5 that 
¥\\li = Sl 5lds = -\l SSss ds 

sc 2 aS\ess \ IVss \ds + 2\l S\est \ \fj,t I ds 

2|||a|||_2,oPIII-i,o||e||2,2||rj||i,2 + 2p|||_i,olkt||i,i||rj||o,i- 



Using the fact that |||cr|||_2.o = Ap]i along with the inequality (A.52), we conclude 
||<5||o,i*;2A[a]||e||2,2||rj||i,2 + 2||e,||i,i|l^,||o,i 

< (1 + A[a])^/E^^/EM ^ (E^mf^^VEM. 



(A.53) 



using Lemma |4.8| and Lemma |5.1| 

We compute the energies of e using the same technique as in Theorem |6.1| we 



try to bound in terms of E„i for to = and to = 1. The lowest one is easy: 



stf + a\es\^) ds 



we have by (A. 51) that 

dEo[e] _ d 
dt ~ dt 

= C[<j]Eo[e] + 2 Ss(ri,,et) + 2 S(f],„et) ds, 
since the boundary term vanishes. Then we can estimate the rest: 
dEo[e] 



dt 



^ C[a]Eo[e] + 2|||^|||o,i||<5||o.i||edlo,o + 2 



-i,o||»7lli,2||et||o,o 



(A.54) 



<: 



C[a]Eo[e] + 2||<5||o,iV^BoM(lil?j|||o,i + \\v\\h2)- 



Now from Lemma 



we have ||?7||i.2 a/ E2 [fj] , while the term |||^|||o,i is a bit 
more difficult (since we have no weighting on the supremum and can't use ( 4.11[ )). 
Instead we use the standard Sobolev inequality (4.9): 

lll^lllo,! ^ ll^llo.l + ll^llo,2 ^ EM, 

using Lemma |4.8| This yields 



(A.55) 



We thus cannot bound EqIe] without also bounding this is not surprising 

since it's hard to even make sense of equation (A. 50 1 without ijst and r]ss both being 
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in . Our next step is to perform the same estimates for i?i[£], at which point 
the estimates do close up. 
We can compute that 



dt 



{^\^st?+^^\ess?) ds ^ C[w]\\e,\\l, + 2C[a]A[a]\\. 

+ 2jg aass(£s,£st')ds + 2'\laSss(f]s,£styds 
+ 4='\l^^AVss,£st')ds + 2'\laS(f]^^^,estyds, 



2,2 



(A.56) 



where again the boundary term vanishes. 

We also want to use a trick to simplify the estimates a bit — the same trick we 
used in deriving (A. 46) — to reduce the derivatives on e to compensate for the high- 
derivative term Sss- Since (ej, tJ^) = 0, we have (^est,Vs} + i^s^Vst) = Oj ^^'^ ^^^^e 
this to write l^(7Sss(fjg, Est) ds = — If^aSssiVst^ ^s} ds. Now (A.56 1 becomes 

i Id + ^^'k-l') ds ^ + A[a])E,[s] + 2lla\a,s\\eA\e.t\ ds 

+ 4 Id ^|<5, 1 177,, \\e,t\ds + 2lla\S\ \rj^^^ \\e,t\ds + 2ll a\S,, \ \rj^^ \\e,\ds. 
Using Corollary |4.5[ we can easily bound all but the last term: 



dt 



Jo 



Sstl^ + a^\£ss\^) ds < C[a]il + A[<j])E,[e] 



2A[a]\\et\\i,i (V^aM^ + V^ll^llo,i) + s\Sss\\vJ\es\ ds. 

(A.57) 



We can bound £'3[ct] ^ E;}[rf]^ using Le mma|5.4[ while (A. 53) bounds ||(5||o,i- 

Finally we deal with the last term of ( A.57[ ) by plugging in (A.50) and using the 
fact that l^ssP + |£ssP = ^(l^s^iP + I^P?2n ill order to eliminate the seemingly 
nonlinear dependence on e. After a series of computations as above, we obtain 



J" 

JO 



s\Sss\\rist\\es\ds<E3[r]fEi[e]. 



Therefore (A.57) becomes 
d 



dt 



\ (alestf + a^lsssf) ds < EM^E^iel 
Jo 



and combining this with ( |A.55 ), we obtain 



dEi[e] 
dt 



< 



E^[r]fEi[e]. 



Now using the inequality s ^ B[(T]a{s) for all s, we bound Ei[e] in terms of -Ei[e]: 

E,[e] (1 + B[a]f E,[e]. 

Using B[a] ^ -S[(Ti] + B[a2] and the bound ( 5.10[ ) for i?[cri] and i?[(T2] in terms of 
^2['7i] and -E2['72] respectively, we ultimately find that 

dEi 



dt 



^ N{t)E,{t), 



(A.58) 



where N{t) is a function depending only on the energies £"3 [771] and -E3[772], which 
are uniformly bounded by assumption. 
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Using Gronwall's inequality, we conclude that if Ei{0) = 0, then Ei{t) = for 
all time. In particular we conclude that a{t, s)\es{t, s)\^ ds = for all t e [0, T], 
so that dss{t, s) = for all t e [0, T] and s e [0, 1]. Since e{t, 1) = 0, we must have 
e{t, s) = for all t and s, whence we conclude r]i{t, s) = r]2it, s) for all t and s. The 
fact that (71 = 0-2 follows. □ 
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